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ABSTRACT 


Research  related  to  the  detection,  location,  and  Iden¬ 
tification  of  underground  nuclear  explosion  devoted  to 
several  major  problem  areas  has  been  completed.  Effects  of 
earth  structure  and  focal  mechanism  on  body-wave  magnitude 
determinations  have  been  investigated,  and  empirical  relations 
between  nuclear  yield  and  Rayleigh  wave  magnitudes  established 
for  selected  seismograph  stations.  Calculations  of  crustal 
transmission  coefficients  for  P-waves  have  been  extended  to 
take  into  account  dip  and  the  results  applied  to  observations 
of  laboratory  explosion  and  deep  earthquakes.  A  model  for 
the  triggering  of  earthquakes  by  explosions  based  on  finite- 
strain  theory  has  been  suggested,  the  superposition  of  the 
explosion-generated  and  triggered-earthquake-generated  P-wave 
computed  for  single  models,  and  aftershocks  of  the  Benham 
event  analyzed.  Methods  of  using  P-wave  and  Rayleigh  wave 
spectra  for  determining  earthquake  mechanisms  have  been  de¬ 
veloped  and  applied.  A  small  local  network  for  detecting 
and  locating  low  magnitude  earthquakes  in  southeast  Missouri 
has  yielded  good  data  on  current  local  seismicity. 
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1.  On  the  Problems  of  Determining  Body-Wave  Magnitude 

Otto  W.  Nuttli 

One  of  the  most  successful  means  of  discriminating  be¬ 
tween  earthquakes  and  underground  explosions  compares  the 
body  and  surface  wave  magnitudes  of  the  event.  This  discrim¬ 
ination  method  was  developed  for  larger  magnitude  events  (mb 
between  about  5  and  6.^',  and  is  gradually  being  extended  to 
smaller  ones. 

At  least  two  classes  of  problems  arise  when  smaller 
magnitude  events  are  considered .  The  first  concerns  the 
physics  of  the  source  mechanism.  For  earthquakes  of  small 
magnitude  the  volume  of  the  source  region  and  the  time  dura¬ 
tion  of  the  source  become  small,  so  that  the  earthquake  source 
will  more  closely  approximate  an  impulse  in  time  applied  at 
a  point.  For  this  reason  the  differences  in  seismic  signals 
from  earthquakes  and  underground  explosions  are  expected  to 
become  less  pronounced  as  their  magnitudes  become  smaller. 

The  second  class  of  problems  peculiar  to  small  magnitude 
events  results  from  the  fact  that  the  amplitudes  of  the  seis¬ 
mic  signals  are  smaller,  so  that  they  are  recorded  by  perhaps 
only  a  few  seismograph  stations  or  arrays.  As  a  consequence 
certain  factors  which  affect  the  amplitude  of  the  recorded 
signal  take  on  a  new  significance.  For  the  larger  magnitude 
events  these  factors  can  to  a  large  extent  be  ignored,  by 
assuming  that  they  randomly  affect  the  ground  motion  at  a 


lurv.c  numtfi*  of  ;;o pj.  stations.  Although  the  scatter 
In  magnitude  values  1’rom  individual  stations  will  be  large, 
the  statistical  average  of  all  such  magnitudes  appears  to 
give  a  reasonable  approx lmat ion  to  the  actual  magnitude. 
However,  when  data  from  only  a  few  stations  are  available, 
the  effects  of  crustal  structure,  azimuthal -dependent  radia¬ 
tion  pattern  of  the  source  and  the  attenuation  of  wave  motion 
with  distance  take  on  much  more  importance. 

In  another  paper  in  this  report  (Effect  of  Regional  Cor¬ 
rection  on  the  Value  of  mb),  Syed  offers  a  methodology  for 
correcting  for  the  azimuthal-dependent  radiation  pattern  of 
the  source.  His  procedure,  which  can  be  applied  in  a  routine 
manner,  does  not  require  that  a  focal  mechanism  solution  be 
determined  for  each  event.  Furthermore,  his  study  enables 
him  to  determine  which  seismograph  stations  will  receive 
anomalously  small  amplitude  P  waves  from  earthquakes  in  a 
particular  geographic  region  because  of  the  focal  mechanism 
pattern  of  the  region.  For  example,  LASA  will  have  smaller 
than  average  amplitude  P  waves  from  earthquakes  beneath  the 
Aleutian  Islands,  but  normal  or  larger  than  average  amplitude 
P  waves  for  earthquakes  beneath  the  Aleutian  Trench. 

The  effects  of  the  crust  and  upper  mantle  structure  on 
the  ground  motion  can,  in  theory,  be  compensated  by  using 
the  crustal  transfer  function.  These  effects  include  the 
frequency-dependent  constructive  and  destructive  interference 
of  waves  reflected  at  the  crustal  boundaries  and/or  converted 
from  P  to  SV  motion.  However,  in  practice,  the  inverse 


problem  of  determining  crustal  structure  from  observations 
of  the  crustal  transfer  function  only  has  been  successful 
for  jng-p?rlod  P  waves.  Attempts  at  using  the  short-period 
P  spectral  ratios  for  crustal  structure  determinations  have 
proved  inconclusive  or  unsatisfactory  (Utsu,  1900,  and  Hase- 
gawa,  1970).  This  can  be  attributed  In  great  part  to  the 
fact  that  the  maxima  and  minima  of  the  crustal  transl’er  func¬ 
tion,  when  considered  a3  a  function  of  earth  period,  are  much 
more  closely  spaced  together  at  periods  of  about  0.*3  to  2.0 
seconds  than  at  periods  of  about  ‘j  to  10  seconds. 

The  extent  to  which  crustal  structure  can  influence  the 

short-period  P  wave  ground  motion  may  be  seen  in  the  plot  of 

the  maximum  A/T  for  the  vertical  component  of  the  P  motion 

produced  by  L0N0SH0T  (Lambert  et  al,  1970,  fig.  4).  At  dis- 
o  o  . 

tances  of  20  to  100  the  spread  in  A/T  values  at  any  given 
distance  is  at  least  one  and  as  much  as  two  orders  of  magni¬ 
tude  (differences  of  10  to  100  times).  As  a  second  example, 
the  data  in  Table  1,  which  contains  mfa  values  from  short- 
period  PZ  for  two  Soviet  underground  explosions  at  Kazakh, 
indicate  a  spread  of  one  magnitude  unit  between  m^  values  at 
the  various  WWSSN  stations.  However,  at  no  station  does  the 
difference  in  mb  for  the  two  events  exceed  0.27  magnitude 
units,  the  average  value  of  the  absolute  difference  being 
only  0.12.  This  indicates  that  the  variation  in  magnitude 
values  between  stations  for  a  single  event  is  not  caused 
principally  by  noise  or  errors  in  reading,  but  rather  is 
reproducible.  Such  reproducibility  is  to  be  expected  if  the 
observed  amplitudes  depend  strongly  upon  the  crustal  structure. 


' fiia t  : i.l  unp.-:*  nun-iv  :H.ru  twiv  him  u  leaser  modify¬ 
ing  ci'iV-'.  upon  ihc  irnplltude  of  th«*  incident  P  wave  for 
the  lonr -period  wjv»m,  aa  -ompared  to  the  short-period  ones. 
PI, '.are  1.1  .'.Ives  the  A/T  values  of  ? for  the  underground 
explosion  at  Novaya  .lernlaya  on  Oct.  iV,  1966.  The  circles 
correspond  to  the  maximum  value  of  A/T  in  the  first  three 
cycles  of  the  P  motion,  as  obtained  from  the  short-period 
WWoSN  seismograms .  "he  crosses  correspond  to  the  A/T  value 
of  the  first  half-cycle  of  the  P  motion  as  recorded  by  the 
long-period  WW33N  seismographs.  The  solid  line  is  the  ex¬ 
pected  value  ol  A/T,  based  on  Gutenberg  and  Richter's  (1956) 
magnitude  :hart,  for  mt  =6.4.  Note  that  in  general  there 
is  more  scatter  amorg  the  short-period  values  than  the  long- 
period  ones.  There  are  three  anomalously  large  values  of 
the  long-period  data  between  65°  and  70°.  These  anomalous 
amplitudes,  which  are  from  stations  on  islands  in  the  western 
Pacific,  will  be  discussed  later. 

Similar  examples  can  be  given  for  LASA,  with  large 
variations  In  short-period  P  wave  amplitude  across  the  array 
(Broome  et  al,  1967,  fig.  A-l)  but  practically  no  variations 
in  long-period  P  wave  amplitudes  across  the  array  (Syed,  1969) 
Ben-Menahem,£t  a_l  (1965,  fig.  9)  gave  the  impulse  re¬ 
sponse  of  a  long-period,  vertical  component  seismograph  system 
(Ts  =  30  sec.,  Tg  *  100  sec.)  to  an  incident  P  wave  for  a  half 
space  model  and  five  models  including  a  layered  crust  over  a 
half-space.  The  crustal  models  include  extreme  types,  such 
as  thin  oceanic  and  thick  mountainous  ones.  For  these  two 


extreme  models  the  P  wave  first  hall  amplitude  was  l.rj  to 
2  times  larger  than  that  of  the  half-space  (O.g  to  0.-;  of  an 
order  of  magnitude)  whereas  for  the  other  three  models  the 
amplitudes  were  almost  the  mine  ar  far  the  haif-space.  This 
figure  not  only  helps  to  exploit;  t h< ■  three  anomalously  large 
amplitudes  of  P  for  the  Novaya  .'.etniya  ev«  at  as  recorded  at 
the  island  stations  at  distances  of  v,°  to  7C°,  but  also 
serves  to  explain  why  the  scatter  of  >/T  values  of  long-period 
P  waves  is  relatively  small. 

Eune,  et  al  (1970)  note  that  m^  values  of  earthquakes  as 
determined  by  the  U.S.C.G.S.  are  consistently  smaller  than 
those  determined  by  the  seismological  stations  of  the  U.S.S.R. 
The  differences  are  sometimes  as  large  as  0.7  magnitude  unit. 
They  believe  the  discrepancies  are  caused  by  differences  in 
the  frequency  response  of  the  seismograph  systems,  the  U.S.C.G. 
values  being  based  upon  data  from  peaked  short-period  systems. 

All  the  reasons  mentioned  above  argue  for  mb  determina¬ 
tions  based  upon  A/T  values  of  the  first  half-cycle  of  the 
long-period  P  wave  motion.  At  present  this  is  only  oractical 
for  earthquakes  or  explosions  with  mb  >  6.0,  because  the  P 
waves  of  smaller  magnitude  events  as  recorded  by  WVv'SSN  seis¬ 
mographs  (most  with  a  peak  magnification  of  1300  or  3000)  are 
too  small  to  be  seen  on  the  seismograms.  However,  the  new 
hlgh-gain,  long-period  seismographs  (Pomeroy  et  al,  1969) 
which  can  be  operated  at  peak  magnifications  of  several  hun¬ 
dred  thousand  should  allow  P  waves  to  be  detected  at  tele- 
seismic  distances  for  mb  *  ^-0- 


It  l.i  :v  ^n:r.«  :«!cJ  that,  whenever  pooalblu,  body  w:r.  •• 
magnitude  v. i  i  .<  n  ic  determined  from  1  rn /.-period  au  we  1 1  ju 
anort-perlcl  P  h.ivc  data.  The  long-period  magnitude  might 
Le  .•ailed  xn|5 .  to  distinguish  It  from  the  short-period  magni¬ 
tude.  A  Bet  of  statistics  should  be  accumulated,  relating 
mB  10  also  comparing  the  standard  deviation  o f  mg  to 

that  of  mb .  It  Is  expected  that  the  nig-M^  relation  will 
not  show  ouch  a  strong  regional  dependence  as  has  been  found 
for  the  mb-M3  relation.  Furthermore,  it  is  expected  that 
the  standard  deviation  of  will  be  less  than  that  of  m^. 

In  magnitude  determinations  the  effect  of  attenuation 
on  the  P  wave  amplitude  appears  in  the  Q  te^m,  where 

mb  .  Q(  A  ,h)  +  log  AA 

and  A  is  eplcentral  distance,  h  focal  depth,  A  the  amplitude 
of  the  ground  motion  and  T  its  period.  All  seismologists, 
except  those  In  eastern  Europe  and  the  Soviet  Union,  presently 
use  the  Q  values  given  by  Gutenberg  and  Richter  (1956)  for 
teleseismic  distances.  For  smaller  distances  Evernden  (1967) 
has  shown  that  lateral  variations  in  crust  and  upper  mantle 
structure  are  large  enough  so  that  different  Q  functions  are 
required  for  different  geographic  regions  at  distances  less 
than  about  2500  Km.  Soviet  seismologists  use  the  Q  function 
of  Vanek  et  ai  (1962)  for  their  mb  determinations.  Recently 
Duda  (1970)  has  presented  a  new  set  of  Q(  A  ,h)  values  for  P 
waves,  based  on  the  assumption  that  the  P  wave  attenuation 
is  caused  only  by  geometrical  spreading  (effects  of  anelas- 
ticity  are  negligible)  and  utilizing  the  1968  P  Travel-Time 
Tables  of  Herrin  et  al  (1968). 
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Pig.  l.T  compares  A/T  values  of  the  vertical  component 
of  P  for  an  event  with  mb  ■  6.4,  using  the  Gutenberg-Richter, 
Vanek  et  al  and  Duda  y  functions,  as  well  as  a  fourth  de¬ 
rived  from  the  Jeffreys-Bullen  travel-time  tables.  All  the 
curves  were  made  to  coincide  at  6o°.  The  greatest  differ¬ 
ences  between  the  curves  occur  from  about  26°  to  45°,  where 
they  can  be  as  large  as  the  equivalent  of  0.6  magnitude  unit. 
Between  45°  and  93°  the  differences  are  smaller,  never  exceed¬ 
ing  0.3  magnitude  unit  and  more  commonly  about  0.1  unit. 

If  sufficient  data  are  available  over  a  large  range  of 
eplcentral  distances,  all  four  calibration  functions  will 
give  approximately  the  same  magnitude  for  a  given  event. 
However,  if  only  a  few  amplitude  data  are  available  and  par¬ 
ticularly  if  they  are  at  distances  less  than  46°,  then  the 
choice  of  the  calibration  function  will  have  a  greater  effect 
on  the  computed  body-wave  magnitude.  The  question  thus 
arises:  Which  of  the  various  magnitude  calibration  functions 

is  closest  to  the  one  that  applies  to  the  real  Earth? 

Gutenberg  and  Richter's  Q  function  is  based  on  theoreti¬ 
cal  calculations  using  their  travel-time  tables,  the  calcu¬ 
lated  values  then  being  modified  to  conform  better  to  actual 
amplitude  data.  There  is  a  very  large  amount  of  scatter 
present  in  the  amplitude  data  they  used  (see,  for  example, 
Gutenberg,  1938,  fig.  1),  so  much  so  that  the  existence  of 
the  majority  of  maxima  and  minima  in  the  curve  seems  ques¬ 
tionable.  The  amplitude  data  also  scatter  widely  about  the 
Q  curve  of  Vanek  et  al  (Bune  et  al,  1970,  fig.  7).  These 


authors,  however,  elected  to  construct  a  more  smooth,  mono- 
tonically  decreasing  amplitude-distance  curve;  an  exception 
to  this  statement  is  a  shallow  minimum  in  the  curve  between 
63°  and  70°.  The  amplitude-distance  curves  calculated  from 
the  smoothed  Jeffreys-Bullen  and  the  Herrin  et  al  tables 
(Duda,  1970)  also  are  generally  monotonically  decreasing 
functions  of  distance;  exceptions  occur  where  there  is  a 
greater  than  average  curvature  of  the  travel-time  curve  (large 
amplitude)  or  a  lesser  than  average  curvature  (small  ampli¬ 
tude  ) . 

Nuclear  explosions,  with  their  almost  symmetrical  P 
wave  radiation  pattern,  have  been  recognized  to  be  a  valu¬ 
able  source  of  data  for  Improving  the  Q  function  or  the 
amplitude-distance  curve.  Carpenter  et  al  (1967)  and  Kalla 
(1970)  constructed  such  curves  using  A/T  values  of  short- 
period  PZ  waves  from  large  underground  explosions.  In  gen¬ 
eral  the  scatter  In  their  data  is  large,  which  can  best  be 
seen  by  comparing  the  two  radically  different  amplitude- 
distance  irves  which  they  derived  as  best  fitting  the  data. 

To  avoid  the  problem  of  scatter  which  appears  to  be 
inherent  when  using  short-period  amplitude  data,  we  have 
measured  A/T  values  of  PZ  recorded  by  long-period  WWSSN 
seismographs  for  underground  explosions  at  Novaya  Zemlya 
(Oct.  27,  1966),  Nevada  (GREELEY,  Dec.  20,  1966  and  BOXCAR, 
Apr.  26,  1968)  and  Amchltka  (MILROW,  Oct  2,  1969) .  Pig.  1.3 
presents  the  data.  In  it  the  MILROW  values  have  been  scaled 
down  by  0,5  magnitude  unit  to  make  them  comparable  to  the 
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values  for  the  other  three  explosions,  all  of  which  are 
taken  to  have  a  body-wave  magnitude  of  6.4.  The  bigger 
A/T  values  of  MILROW  are  caused  not  only  by  larger  ampli¬ 
tudes,  but  also  by  shorter  periods.  Crosses,  circles  and 
triangles  serve  to  distinguish  data  from  Novaya  Zemlya, 
Nevada  and  Amchltka,  respectively.  Small  symbols  Indicate 
less  reliable  values,  resulting  either  from  the  presence 
of  microseismic  noise  or  the  fact  that  the  trace  amplitude 
of  PZ  was  less  than  0.8  mm.  Five  data  points  with  anoma¬ 
lously  large  A/T  values  have  been  omitted  from  the  figure. 
Two  curves  are  included,  the  one  the  Gutenberg-Richter  for 
m^  =  6.4,  the  other  a  curve  drawn  more  or  less  intuitively 
as  a  "by  eye"  fit  to  the  data.  Differences  between  the 
two  curves  are  greatest  from  40°  to  70°,  where  they  are 
as  large  as  0.4  magnitude  unit. 

Both  amplitude  curves  of  Fig.  1.3  are  empirical  ones 
which  have  been  constructed  in  a  rather  arbitrary  manner. 

To  this  point  in  the  discussion  they  have  not  yet  been 
subjected  to  the  constraints  Imposed  upon  them  by  the 
travel-time  curve. 

For  an  event  with  zero  focal  depth,  the  amplitude  A 
beneath  the  station  is 

A  *=  C  ^|di/d&  Jtan  i/sin~A_' 

where  C  is  a  constant,  i  is  the  angle  of  Incidence  at  the 
focus  and  the  epicentral  distance.  The  earth  model  is 
assumed  to  be  spherically  symmetric  and  the  effects,  if  any, 
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of  anelasticity  are  ignored.  The  amplitude  Az  of  the  verti¬ 
cal  component  of  the  surface  P  motion  is  related  to  A  by 

Az  =  A  cos  T 

where  i  is  the  apparent  angle  of  incidence  at  the  station. 

Nuttll  (1964)  experimentally  determined  the  relation  of  i 
to  A  for  P  waves  of  various  periods . 

The  slope  of  the  travel -time  curve,  dT/dA  >  for  a  sur¬ 
face  focus  is 

dT/d A  =  r  sin  i/v 

where  r  is  the  radius  of  the  earth  model,  v  the  velocity  and 
i  the  angle  of  incidence  at  the  surface. 

The  three  equations  above  provide  implicit  relations 
between  the  dependent  variables  Tz,  i,  dT/dA  and  the  inde¬ 
pendent  variable  A  .  Because  the  dT/dA  vs  A  relation  for  P 
waves  is  known  with  a  high  degree  of  accuracy,  the  i  vs  A  and 
vs  A  relations  are  severely  constrained.  For  example,  it 
was  found  by  calculations  based  on  the  above  equations  that 
the  amplitudes  of  both  empirical  curves  of  Fig.  1.3  are  too 
large  between  about  35°  and  45°.  As  drawn  they  would  imply 
errors  of  up  to  five  or  more  seconds  in  the  travel-time  curve 
at  those  distances,  which  is  obviously  not  acceptable. 

Detailed  calculations  show  that  iruch  of  the  "topography," 
that  is,  the  existence  of  pronounced  maxima  and  minima  in 
the  amplitude-distance  curves  of  Fig.  1,3>  is  not  compatible 
with  the  Jeffreys-Bullen  or  the  Herrin  et  al  P  tables.  From 
this  it  follows  that  the  Gutenberg -Richter  magnitude  calibration 
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function  for  zero  focal  depth  is  inconsistent  with  observed 
travel  times. 

A  trial  and  error  method  was  used  to  calculate  an  ampli¬ 
tude-distance  curve  which  satisfied  the  travel-time  tables 
(maximum  permissible  difference  of  one  second)  and,  as  well 
as  possible,  the  amplitude  data  of  Pig.  1.3.  Prom  this  curve 
a  Q  magnitude  calibration  function  was  calculated  and  tabulated. 
Its  values  are  given  in  Table  2. 

Pig.  1.4  compares  the  new  Q  function  with  those  of  Guten¬ 
berg  and  Richter  and  of  Duda .  In  general  it  agrees  well  with 
Duda's,  their  separation  being  everywhere  less  than  or  equal 
to  0.1  magnitude  unit  except  from  27°  to  28°,  where  it  is 
about  0.2  unit.  Differences  between  the  new  Q  function  and 
that  of  Gutenberg  and  Richter  are  larger,  being  as  big  as 
0-2  magnitude  unit  from  37°  to  44°,  47°  to  48°  and  62°  to 
69  .  The  general  shape  of  the  new  Q  curve,  or  the  related 
A/T  curve,  conforms  with  Duda's  assumption  that  the  effects 
of  anelastlclty  are  small  compared  to  those  of  geometrical 
spreading. 

Table  3  gives  magnitudes  of  four  underground  explosions, 
using  the  amplitudes  of  PZ  recorded  by  long-period  WWSSN 
seismographs.  Magnitudes  were  calculated  using  the  new  Q 
function  and  that  of  Gutenberg  and  Richter.  Por  three  of 
the  four  events  the  body-wave  magnitude  calculated  using  the 
new  Q  function  has  a  slightly  smaller  standard  deviation. 

Por  the  fourth  event,  BOXCAR,  the  standard  deviations  are 


the  same . 
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Table  4  compares  magnitudes  of  two  Aleutian  earthquakes, 
using  amplitude  data  of  Syed  (personal  communication)  which 
have  been  corrected  for  focal  mechanism.  Average  magnitudes 
and  their  standard  deviations  are  given,  for  magnitudes  cal¬ 
culated  using  the  new  Q  function  as  well  as  those  of  Duda  and 
of  Gutenberg  and  Richter.  Prom  the  table  it  can  be  seen 
that  the  average  magnitudes  and  standard  deviations  obtained 
using  the  Duda  and  new  Q  functions  are  the  same,  and  they 
differ  only  slightly  from  those  based  on  the  Outenberg-Richter 
Q  function. 


Summary 

The  purposes  of  this  paper  are  to  examine  the  factors 
which  cause  problems  in  body-wave  magnitude  determination,  to 
assess  their  relative  importance  and,  if  possible,  to  suggest 
means  of  reducing  the  errors  in  such  evaluations. 

Three  principal  factors  which  can  cause  such  problems 
are  the  azimuth-dependent  radiation  pattern  of  the  earthquake 
source,  the  effects  of  crust  and  upper  mantle  structure  at 
the  station  and  an  inexact  knowledge  of  the  magnitude  cali¬ 
bration  function. 

The  new  magnitude  calibration  function  which  has  been 
developed  in  this  study  differs  so  little  from  that  of  Duda 
(1970)  that  the  two  may  be  considered  to  be  the  same.  They 
are  believed  to  be  an  improvement  over  the  Outenberg-Richter 
(1956)  function  because  they  are  more  consistent  with  observed 
travel  times.  The  latter  function  will  underestimate  m^  by 
about  0.2  magnitude  unit  between  37°  and  44°  and  overestimate 
it  by  about  the  same  amount  between  62°  and  69°.  Otherwise 
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the  three  calibration  functions  will  give  almost  the  same 
value  of  mb. 

Ignoring  the  effects  of  the  focal  mechanism  of  earth¬ 
quakes  can  lead  to  discrepancies  as  much  as  0.5  magnitude 
unit  at  particular  stations,  depending  upon  the  radiation 
pattern.  However,  Syed  has  developed  a  straightforward 
method  for  taking  account  of  this  phenomenon  which  is  readily 
adaptable  to  routine  calculations,  so  that  this  problem  can 
be  said  to  be  nearly  solved. 

The  biggest  source  of  error  in  mb  determination  is  the 
effect  of  the  crust  and  upper  mantle  structure  at  the  station. 
For  short-period  P  waves  it  can  cause  errors  of  as  large  as 
one  magnitude  unit  at  particular  stations.  Unfortunately  the 
station  correction  is  a  complicated  function  of  epicentral 
distance,  the  frequency  content  of  the  incident  wave,  the 
geological  structure  at  the  station  and  the  azimuth  of  wave 
approach  if  there  is  lateral  heterogeneity  in  geological 
structure  near  the  station.  Probably  the  only  way  of  arriv¬ 
ing  at  station  corrections  for  short-period  data  is  on  a 
statistical  basis,  which  is  a  formidable  task.  However,  the 
use  of  long-period  P  data  would  avoid  many  of  these  diffi¬ 
culties.  For  such  data  the  largest  error  to  be  expected  at 
a  particular  station  because  of  crustal  structure  is  only 
about  0.3  magnitude  unit.  Furthermore,  because  the  effect 
of  the  crust  and  upper  mantle  structure  on  the  amplitude  of 
long-period  P  waves  requires  but  a  coarse  knowledge  of  the 
structure,  simple  calculations  can  be  made  to  estimate  this 
effect  at  individual  stations. 
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By  using  the  first  half-amplitude  of  long-period  P,  after 
correcting  it  for  crustal  structure  and  focal  mechanism, 
better  estimates  of  body-wave  magnitude  can  be  expected.  It 
seems  likely  that  the  standard  deviation  of  m^  can  be  made  as 
small  as  0.2  and  possibly  0.1  magnitude  unit. 
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STA. 

A 

0 

Jan. 15,  1965 

March  20,  1966 

AAM 

86.9 

6.13 

5.86 

0.27 

ALQ 

95.7 

5.67 

5.64 

0.03 

ATL 

95.7 

5.96 

5.88 

0.08 

BEC 

92.0 

- 

5.77 

- 

BKS 

9C.9 

5.85 

6.06 

0.21 

BLA 

91.3 

5.71 

5.77 

0.06 

BOZ 

84.6 

5.99 

6.17 

0.18 

CMC 

62.3 

6.04 

6.19 

0.15 

COL 

60.1 

6.52 

6.33 

0.14 

COR 

84.3 

6.51 

DAL 

97.6 

6.33 

- 

- 

DUG 

90.3 

- 

5.92 

- 

FLO 

91.5 

5.72 

5.57 

0.15 

ODH 

56.0 

- 

5.98 

- 

GEO 

89.0 

6.02 

5.76 

0.26 

GOL 

91.0 

5.87 

5.89 

0.02 

osc 

94.7 

6.02 

5.88 

0.14 

JCT 

100.1 

- 

5.59 

- 

LON 

82.4 

- 

6.07 

- 

LUB 

97.1 

- 

6.16 

- 

MDS 

86.6 

6.20 

6.22 

0.02 

MNN 

85.4 

5.98 

- 

- 

OGD 

86.4 

5.71 

- 

- 

OXF 

95.6 

6.20 

6.23 

0.03 

RCD 

86.7 

5.98 

6.09 

0.11 

SCP 

87.3 

- 

5.81 

TUC 

98.0 

5.59 

5.48 

0.11 

Average 

6.00  +  0.26 

5.93  +  0.24 

0.12 

Table  2.  Q  Values  for  PZ  (Zero  Focal  Depth) 
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A 

Q 

A 

Q 

A 

Q 

20 

6.21 

48 

6.67 

76 

6.86 

21 

6..  *3 

49 

6.69 

77 

6.86 

22 

6.26 

50 

6.71 

78 

6.86 

23 

6.28 

51 

6.75 

79 

6.86 

24 

6.32 

52 

6.78 

80 

6.87 

23 

6.37 

53 

6.80 

81 

6.87 

26 

6.44 

54 

6.82 

82 

6.88 

2? 

6.63 

55 

6.79 

83 

6.89 

28 

6.76 

56 

6.76 

84 

6.90 

29 

6.78 

57 

6.74 

85 

6.91 

30 

6.77 

38 

6.72 

86 

6.92 

31 

6.72 

59 

6.73 

87 

6.92 

32 

6.71 

60 

6.75 

88 

6.93 

33 

6.71 

6l 

6.76 

89 

6.94 

34 

6.71 

62 

6.78 

90 

6.97 

35 

6.71 

63 

6.79 

91 

7.01 

36 

6.71 

64 

6.81 

92 

7.06 

37 

6.63 

65 

6.83 

93 

7.11 

38 

6.65 

66 

6.85 

94 

7.15 

39 

6.65 

67 

6.85 

95 

7.20 

40 

6.65 

68 

6.83 

96 

7.24 

41 

6.65 

69 

6.81 

97 

7.29 

42 

6.65 

70 

6.80 

98 

7.34 

43 

6.65 

71 

6.81 

99 

7.38 

44 

6.65 

72 

6.83 

100 

7.43 

45 

6.65 

73 

6.85 

46 

6.65 

74 

6.86 

47 

6.65 

75 

6.86 
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Table  3.  rnB  from  Long-Period  PZ  (Explosions) 


Location 

Date 

New  Q 

mB 

Gut-Rich. Q 

No. 

of 

Sta. 

Novaya 

Zemlya 

Oct.  27,1966 

6.44+0.19 

6.49+0.23 

46 

Nevada 

(GREELEY) 

Dec.  20,1966 

6 . 36+0 . 19 

6.38+0.23 

29 

Nevada 

(BOXCAR) 

Apr.  28,1968 

6.32+0.26 

6.35+0.26 

28 

Amchitka ( MILROW) 

Oct.  2,1969 

6.87+0.20 

6.94+0.27 

14 

20 
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A  (deg.) 


7.5 


A  (deg ) 


2.  Effect  of  Regional  Correction  on  the  Value  of  mb. 

Atlq  Syed 

Discrepancies  In  magnitude  determination  are  known  to 
result  in  part  from  the  effect  of  the  earthquake  radiation 
pattern  on  the  recorded  wave  amplitudes.  The  purpose  of  the 
study  reported  here  is  to  obtain  a  quantitative  estimate  of 
the  effect  of  the  radiation  pattern  on  P  wave  amplitudes  and 
to  develop  a  methodology  for  improving  body  wave  magnitude 
determinations.  The  results  of  the  study  may  also  be  used 
as  an  auxiliary  discriminant  between  earthquakes  and  explo¬ 
sions  . 


Correction  for  Mechanism 

The  P  wave  displacement  at  the  surface  of  a  focal  sphere 
of  radius  R  for  a  double  couple  point  source  is  given  by  the 
expression  . 


u„  - 


(1) 


If  for  a  given  earthquake  the  observed  P  wave  amplitudes 
are  divided  by  the  product  2xy  the  P  amplitudes  will  be 
corrected  relative  to  the  maximum  value  of  the  P  wave  radia¬ 
tion.  For  magnitude  calculations  in  keeping  with  the  assump¬ 
tions  made  in  the  definition  of  magnitude,  a  further  normal¬ 
ization  is  required,  a  reduction  of  the  P  wave  amplitude  to 
the  average  value  of  the  displacement  of  P  over  the  surface 


of  the  focal  sphere, 
from  the  relation 


This  average  amplitude  la  obtained 
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/zyyJ IS 

u,  -  -Jit—  = 

uf  -  OA Z4(up)m#x 

Thus,  the  total  reduction  of  P  wave  amplitudes  to  take  Into 
account  Includes  first  the  application  of  the  2xy  factor  and 
second  multiplication  by  0.424  to  equate  the  P  wave  radiation 
from  a  double  couple  to  that  from  a  spherically  symmetric 
source  of  the  same  strength. 

It  has  long  been  observed  that  earthquake  foci  located 
In  a  given  region  tend  to  have  similar  orientations  of  the 
focal  mechanism.  This  observation  la  In  keeping  with  the 
theory  of  plate  tectonics.  Plate  tectonics.  In  turn,  provides 
a  basis  for  prediction  of  the  dominant  focal  mechanism.  On 
the  hypothesis  that  characteristic  foci  do  exist  for  given 
hypocentral  regions,  P  wave  amplitudes  can  be  corrected  for 
the  source  mechanism  and  thereby  magnitude  estimates  can  be 
obtained  which  have  smaller  standard  errors. 

A  procedure  for  applying  the  source  mechanism  correction 
la  aa  follows: 

1.  Given  the  source  parameters  for  a  typical  earthquake 
In  a  particular  hypocentral  region,  compute  the  2xy  factor 
for  all  seismograph  stations. 

2.  On  the  basis  of  this  factor,  predict  which  stations 
will  record  relatively  large  P  wave  amplitudes  for  earthquakes 
for  the  region  in  question,  and  which  relatively  small. 


.  7 


3.  Selec*.  z.  *<.  ‘  .  .ne  °  -  r...,niitudes  are 

best  suited  f  c  .•  -urvnquakes  from 

that  region.  .  ■  .  cion  cl  stations 

might  be 

a)  3ta7. .  ..  .  ..  .  t..  *.:.e  ..icentral  distance 

.....  .  .  o  \ r* 

b)  Th-  ..  .  .1  be  greater  than  the 

*  J  4  13  1  1  J.’ j  0  .  "tL'n  . 

The  latter  criterion  insures  that,  the  P  wave  amplitude  be 
relatively  large,  but  also  that  small  differences  between  the 
actual  source  mechanl arc  t/.e  average  «r  cr.aracteristlc 
mechanism  of  the  region  will  ha  'e  -ri.- a. "I'ecc  on  the  re¬ 
duced  amplitudes. 

4  He...-  r.  1:  1  y  t:.  ?  v  .  *•./  _arger  than 

the  norm  a  .id  .  '  i.  i.  the  magni¬ 
tude  determine c  „^.ns  will  be  too  large.  It  is  possible  to 

correct  for  tne  fleet  by  su..  ora  :  t.  lug  from  the  Initial  m^ 

value  the  :v.u:.c.  log  (jxy/ . 424 } ,  -here  2  ~  ;  Is  the  aver¬ 

age  value  of  ..  g  .-11  stations  with  2xy  greater  than  0.424. 
There  is  a  constant,  correction  to  be  applied  to  all  earth¬ 
quakes  In  any  particular  geographic  group.  It  is  called  the 
regional  correction  factor.  In  practice  values  of  this 
correction  in  the  regie .03  studied  thus  far  are  found  to  vary 
from  0.1  to  0.2n  of  a  magnitude  unit. 

5.  If  ittlrablf ,  instead  of  the  required  correction  an 
individual  mee. nanism  correction  may  be  applied  by  subtracting 
the  quantity  log  (2::y/.424)  determined  at  each  station.  This 


correction  is  applicable  even  to  stations  at  which  2xy  -<£.424, 
provided  the  2xy  factors  be  not  too  small.  That  is,  stations 
very  near  nodal  lines  should  always  be  excluded. 

Application 

In  order  to  test  this  suggestion  and  to  exemplify  the 
order  of  magnitude  of  the  corrections  involved,  earthquakes 
from  two  regions  have  been  studied:  earthquakes  of  the 
Aleutian  Islands,  and  ones  of  the  North  Atlantic.  Earth¬ 
quakes  of  the  Aleutian  Islands  (Figure  21)  divide  into  three 
groups  on  the  basis  of  their  source  mechanism  as  determined 
by  Stauder  (1968a, b).  Similarly  those  of  the  Atlantic  divide 
into  two  groups  (Sykes,  1967).  For  each  group  we  can  calcu¬ 
late  the  average  or  dominant  source  mechanism,  and  then  apply 
the  various  corrections  Indicated  above. 

1.  Earthquake  of  March  30,  1965. 

This  earthquake  belongs  to  group  3b  of  the  Aleutian 
Island  earthquakes.  This  group  consists  of  earthquakes 
which  occur  along  a  narrow  line  immediately  below  the  axis 
of  the  Aleutian  trench  or  under  the  seaward  slope  of  the 
trench.  Foci  of  this  group  are  uniformly  extensional  in 
character,  with  the  axis  of  tension  aligned  normal  to  the 
local  axis  of  the  trench.  The  average  focal  mechanism  for 
earthquakes  in  this  group  is  characterized  by  two  nodal 
planes  dipping  about  40°-50°. 

For  n>b  determinations  the  amplitudes  and  periods  of 
the  first  half -cycle  of  the  vertical  component  long  period 


P  waves  were  the  basic  data.  The  2::,r  :  are  determined 

and  listed  in  Table  1  alorv  with  lh  .  :c  d- ter’Vi.ned  the 

ordinary  way  (m^,  unc career.  f.  -  •  t  .at  for  small 

2xy  factors  the  mD  valuer  u:  :  .  tr.jn  for  the  larger 

2xy  factors.  The  averags  rre :-.-.r.isude  using  the 
data  of  all  37  station.-,  *.s  'T.l?  -  l.-.i.  Trie  average  mb  value 
for  all  stations  wi  2.<y  1  r.  -.  I  7-3  n  0.15.  The  aver¬ 
age  value  for  stations  with  2xy  <  0  -+24.  is  7.03  +  0.17. 
Applying  the  regional  focal  mechanism  correction  for  all 
stations  with  2xy  >  0.424,  =  7-12  t  0.15.  Finally,  cor¬ 

recting  for  individual  station  focal  mechmlsm  corrections 
we  obtain  for  all  stations  m^  «=  7.13  j  C.lh.  These  results 
are  summarized  in  Table  2. 

It  is  noteworthy  that  using  oh.  tiers  with  2xy  >  0.24 
and  apply.!*  r,  the  constant  r*3.;io  .  r  section,  the  mb  value 

is  very  close  to  that  "rom  oil  the  stations  using  the  indi¬ 
vidual  mechanism  corrections.  This  v-scTd  seem  to  Justify 
using  for  magnitude  determinations  only  those  stations  for 
which  2xy  >  >  .  One  may,  in  fact,  in  this  way  define 

zones  on  the  earth's  surface,  identifying  stations  (solid 
circles  ir.  figure  2.2)  favorable  for  mb  determinations,  and 
others  (crosses  in  the  figure)  to  be  avoided  in  magnitude 
determinations . 

2.  Earthquake  of  November  22,  1965. 

Thi3  earthquake  belongs  to  group  1  of  the  Aleutian 
Island  earthquakes  .  ThlB  group  consists  of  foci  on  the 
concave  side  of  the  island  arc,  in  the  zone  immediately 


south  of  and  under  the  island  chain.  Focal  mechanism  sta¬ 
tions  in  this  group  have  one  steeply  dipping  nodal  plane, 
with  the  other  nodal  plane  nearly  horizontal.  Many  stations 
in  this  case  are  close  to  the  nodal  plane.  Consequently  the 
number  of  stations  suitable  for  P  wave  magnitude  determina¬ 
tion  was  much  smaller  than  in  group  3b.  The  average  uncor¬ 
rected  magnitude  using  the  data  of  all  19  stations  is 
6.20  +  0.24.  The  average  mb  for  all  stations  with  xy  larger 
than  0.424  corrected  for  the  regional  effect  is  5-99  +  0.22. 
The  average  mb  of  all  stations  whose  amplitudes  have  been 
corrected  for  the  source  mechanism  is  6.00  +  0.26. 

3.  Earthquake  of  August  3*  1963. 

This  earthquake  belongs  to  group  1  of  the  North 
Atlantic  earthquakes.  The  events  in  this  group  are  located 
on  the  equatorial  fracture  zone  of  the  mid-Atlantic  ridge. 
The  mechanism  of  these  earthquakes  correspond  to  two  near 
vertical  nodal  planes.  Consequently,  the  2xy  factors  calcu¬ 
lated  for  a  typical  earthquake  in  this  group  at  all  the  sta¬ 
tions  at  a  distance  of  25-100  degrees  are  smaller  than  0.424, 
which  is  the  normalizing  factor.  Consequently  it  is  not 
possible  to  satisfy  the  second  criterion  in  selecting  the 
stations  suitable  for  P  wave  magnitude  determination.  Depart 
ing  from  this  criterion,  we  have  selected  those  stations 
whose  2xy  factors  calculated  for  the  average  focal  mechanism 
solution  are  larger  than  0.2  and  established  the  appropriate 
regional  corrections.  Since  the  uncorrected  P  wave  magni¬ 
tudes  using  the  data  from  these  stations  will  be  too  low, 


•  1 

the  regional  correction  has  to  be  added  to  the  average  mag¬ 
nitude  value.  The  results  of  the  m^  calculation  for  this 
earthquake  show  that  the  average  uncorrected  magnitude  using 
the  data  of  all  31  stations  is  6.51  +  0.29.  The  average  mb 
for  15  stations  whose  2xy  factors  are  larger  than  0.2  cor¬ 
rected  for  the  regional  effect  is  6.67  +  0.20,  the  average 
m^  for  the  same  stations  whose  amplitudes  have  been  corrected 
for  the  source  mechanism  is  6.73  +  0.18. 

Conclusions 

Prom  this  study  It  can  be  concluded  that  a  significant 
part  of  the  scatter  in  the  determination  of  earthquake  magni¬ 
tudes  is  due  to  the  effect  of  the  earthquake  radiation 
pattern  on  the  recorded  P  wave  amplitudes.  The  correction 
for  the  focal  mechanism  becomes  more  important  when  the 
focal  mechanism  solution  is  one  in  which  one  or  both  nodal 
planes  are  steeply  dipping.  In  these  cases  the  average  mag¬ 
nitude  determined  by  using  the  uncorrected  P  wave  amplitude 
data  may  be  in  error  by  as  much  as  a  quarter  of  a  magnitude 
unit. 

Since  the  observational  data  here  examined  have  shown 
that  the  average  mb  for  stations  with  smaller  2xy  factors 
is  significantly  smaller  than  the  average  m^  for  stations 
with  larger  2xy  factors,  there  exists  a  possibility  of  com¬ 
paring  the  values  of  uncorrected  m^  from  these  two  sets  of 
stations  in  a  region  where  the  average  focal  mechanism 
solution  is  known  to  differentiate  between  the  earthquakes 


and  explosions.  The  radiation  pattern  due  to  the  explosive 
source  being  azlmuthally  uniform,  the  tvo  sets  of  values 
should  be  identical  for  explosion  and  different  for  the 

earthquake . 
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Table  1.  Earthquake  of  March  30,  1965. 


Sta . 

Dlat. 

2xy 

^(uncorr) 

iD^(corr 

LON 

39.0 

0.250 

6.88 

7.11 

CMC 

36.2 

0.256 

6.89 

7.11 

KIP 

3^.7 

0.278 

7.03 

7.21 

DUO 

48.1 

0.322 

7.02 

7.14 

QUA 

45.6 

0.322 

6.79 

6.91 

RCD 

51.5 

0.338 

7.08 

7.18 

OSC 

49.8 

0.346 

7.06 

7.15 

GOL 

52.9 

0.352 

7.14 

7.22 

RAB 

58.9 

0.356 

7.09 

7.16 

HNR 

61.8 

O.360 

7.30 

7.37 

TUC 

54.4 

0.382 

6.93 

6.97 

LUB 

59.0 

0.422 

6.88 

6.88 

MDS 

59-4 

0.424 

7.37 

7.37 

FLO 

62.2 

0.448 

7.34 

7.32 

A  AM 

63.4 

0.466 

7.25 

7.21 

SHA 

69.3 

0.506 

7.57 

7.49 

WES 

69.5 

0.532 

7.05 

6.95 

RIV 

87.3 

0.594 

7.50 

7.35 

TAU 

96.9 

0.624 

7.23 

7.06 

KTG 

58.4 

O.636 

7.57 

7.39 

ADE 

91.9 

O.638 

7.42 

7.25 

BHP 

91.0 

0.642 

7.35 

7.17  . 

QUI 

9  8.8 

0.646 

7.17 

6.99 

Table  1.  (cont'd) 
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§ta  a 

Dlst . 

2xy 

(uncorr) 

/nb(  corr) 

CAR 

97.3 

0 , 664 

7.20 

7.00 

TRN 

100.3 

0.672 

7.13 

6.93 

VAL 

77.6 

0.674 

7.33 

7.07 

MUN 

98.6 

C.714 

7.40 

7.18 

HKC 

56.5 

0.744 

7.25 

7.01 

KEV 

58.1 

0.780 

7.43 

7.16 

KON 

69.7 

0.804 

7.26 

6.99 

PTO 

88.5 

0.826 

7.26 

6.97 

COP 

73.4 

0.834 

7.28 

'  6.99 

TOL 

89-9 

0.844 

7.52 

7.22 

MAL 

93.0 

0.850 

7.17 

6.87 

STU 

80.5 

0.656 

7.26 

6.96 

KOD 

88.7 

0.952 

7.61 

7.26 

JER 

91.6 

O.966 

7.51 

7.15 

<■) 


Table  2.  Average  mb  values,  March  30,  1965. 

Uncorrected 

All  stations  (n  -  37)  mb  *  7*34  ±  0.21 

For  2xy  <.  0.42  (n  *  13)  mb  “  7.03  +  0.17 

For  2xy  >  0.42  (n  -  24)  mb  -  7.3^  +  0.15 

Corrected 

Regional  corr.  (n  -  24)  -  7.12  +  0.15 

Individual  mech.  corr.  (n  -  37)  mb  -  7*13  +  0.16 
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Table  3*  August  3.  1963  Earthquake 


Sta . 

Dlst . 

>r.r 
— - 

•T!-0  v'uncorr ) 

,mb(corr) 

NAI 

36.3 

C.030 

6.42 

7.56 

BOO 

38.2 

C.C36 

5.30 

6.57 

KON 

62.1 

C.GoC 

6.76 

7.60 

NUR 

69.1 

0.062 

6.73 

7.56 

BUL 

69.1 

0 . 0  68 

6.4l 

7.20 

VAL 

48.8 

0.072 

6.73 

7.50 

AAE 

73.7 

0.074 

6.75 

7.51 

PRE 

70.3 

0.084 

6.52 

7.22 

CAR 

30.9 

0.104 

6.05 

6.66 

QUI 

43.3 

0.111 

6.4X 

6.99 

WIN 

59-9 

0.124 

6 . 80 

7.33 

TRN 

23.4 

0.128 

6.38 

6.90 

1ST 

66.0 

0.154 

6.75 

7.19 

ATU 

61.4 

0.180 

6.93 

7.30 

BKS 

82.5 

0.194 

6.17 

6.51 

OSC 

78.4 

0.198 

6.32 

6.66 

TUC 

73.3 

0.212 

6.33 

6.63 

LON 

81.7 

0.214 

6.16 

6.46 

TOL 

42.8 

0.226 

6.65 

6.92 

DUO 

75.2 

0.234 

6.34 

6.60 

ALQ 

69.8 

0.242 

6.53 

6.77 

MAL 

40.6 

0.252 

6 . 70 

6.92 

OOL 

69.4 

0.260 

6.45 

6.66 

Table  3.  (Cont'd) 


Sta . 

Dlst 

2xy 

mb(uncorr) 

®b(corr) 

LPA 

47.3 

0.266 

6.97 

7.17 

NNA 

45.3 

0.268 

6.61 

6.81 

RCD 

68.5 

0.276 

6.39 

6.58 

SHA 

53.9 

0.312 

6.70 

6.83 

FLO 

57.9 

0.322 

6.45 

6.55 

A  AN 

54.4 

0.366 

6.61 

6.57 

OEO 

48.4 

0.398 

6.67 

6.69 

SCP 

49.8 

0.398 

6.69  • 

6.72 

Table  4.  November  22,  1965  Earthquake 


Sta  • 

PMG 

CTA 

MAN 

BAG 

HKC 

TOL 

CHO 

KOD 

HOW 

POO 

KON 

SHL 

JER 

1ST 

SHI 

NDI 

QUE 

TAB 

LAH 


3  b 


Diet. 

2xy 

njk(uncorr) 

mb(corr) 

66.8 

0.124 

5.94 

6.47 

77.0 

0.130 

5.83 

6.34 

59.2 

0.608 

6.50 

6.35 

58.2 

0.628 

6.23 

6.06 

57.7 

0.730 

6.15 

5.92 

89.0 

0.788 

5.85 

5.58 

69.8 

0.798 

6.00 

5.72 

89.6 

0.816 

6.51  ' 

6.22 

73.5 

O.858 

6.38 

6.08 

85.2 

0.860 

6.45 

6.15 

69.0 

0.862 

5.98 

5.67 

69.1 

0.862 

6.11 

5.80 

91.5 

0.886 

6.14 

5.82 

84.5 

O.890 

6.05 

5.73 

86.9 

0.904 

6.49 

6.16 

75.3 

0.904 

6.28 

5.95 

79.5 

0.914 

6.16 

5.83 

82.0 

0.916 

6.17 

5.84 

74.7 

0,920 

6.63 

6.30 
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Earthquake  of  November  22,  1965 

Av.  mb  (unco rr)  -6.20+0.24  n  -  19 

Av.  mfe  (uncorr)  for  2xy  >  0.42  -  6.24  +  0.22  n  -  17 

Av.  mb  (corr  for  regional  effect)for 

2xy  >0.42  -  5.99  +  0.22  r  -  17 

Av.  mb  (corr  for  focal  mech)  -  6.00  +  0.26  n  -  19 

Earthquake  of  August  3,  1963 

Av.  mb  (uncorr)  -  6.51  +  0.29  n  -  31 

Av.  mb  (uncorr)  for  2xy  >  0.20  -  6.55  +  0.20  n  -  15 

Av.  mb  (corr  for  regional  effect)  for 

2xy  >  0.20  -  6.67  +  0.20  n  -  15 

Av.  mfa  (corr  for  focal  mech)  for 

2xy  >0.20  -  6.73  +  0.18  n  -  15 


40 


Figure  1. 


Figure  2. 


Figure  3. 


Titles  to  Figures 

Earthquakes  of  Aleutian  Island  region  divided 
Into  three  groups:  1)  earthquakes  under  the 
islands,  western  portion  of  arc,  2)  earth¬ 
quakes  under  the  islands,  eastern  portion  of 
arc,  3a>b)  earthquakes  with  foci  under  the 
trench. 

Map  indicating  stations  unfavorable  (x)  and 
favorable  (•)  for  determination  of  for 
Aleutian  Island  earthquakes  of  group  3a. 

Map  indicating  stations  unfavorable  (x)  and 
favorable  (f )  for  determining  m^  for  Aleutian 
Island  earthquakes  of  group  1. 
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Nu-l».ar  Yields  from  Rayleigh  Wave3 


Donald  E.  Wagner 

Introduction 

;jurfa  •*?  waves  are  a  function  of  the  depth  of  burst, 
the  time  history  of  the  event,  and  the  wave  guide  through 
whl'h  they  propagate.  For  nuclear  explosions  the  depths 
vary  only  between  l-l  kilometer  while  the  time  history  is 
assumed  constant.  Therefore,  the  wave  guide  is  left  as 
the  factor1  controlling  the  observed  Rayleigh  wave  signature 
at  a  given  station.  Furthermore,  unlike  body  waves,  which 
are  affected  cor.s iderably  by  the  source  medium,  Rayleigh  waves 
are  a  function  of  the  average  velocity  and  density  In  a  layer 
or  wave  guide.  Thus,  they  are  much  less  susceptible  to  in- 
homogenelties  or  velocity  contrasts  both  at  the  shot  point 
and  along  the  path.  The  distinct  similarity  of  nuclear  events 
at  a  station  (see  for  example  Fig.  2,  Shurbet,  1969)  en¬ 
couraged  the  autnor  to  attempt  a  prediction  of  yield  based 
on  Rayleigh  wave  amplitudes. 

In  a  recent  release  (Higgins,  1970),  several  nuclear 
explosions  were  declassified,  offering  a  list  of  known-yield 
events  observable  at  continental  WWSSN  stations.  The  follow¬ 
ing  study  applies  a  revised  (Nuttli,  1970)  surface  wave 
magnitude  formula  based  on  maximum  ratios  of  A/T  from  Ray¬ 
leigh  waves  at  all  stations  to  develop  a  yield  versus  M_ 

s 

(average)  master  curve  for  seven  known  energies.  In  a 


further  attempt  to  reduce  dissimilarities  owing  to  different 
wave  guides,  maximum  ;./?  .  i;;s  were  plotted  as  a  1‘unction 

of  yield  at  eig.nt  s;..::.'i'  . tlur.s.  Ail  in  all,  nine  master 

curves  were  plotter,  from  w.-.ic.g  eleven  unknown  yields  were 
predicted . 

I’ro.  edure 

The  Richter  •...■’face  wave  I'orr.uJ.i  V. „  -  j .  1  +  log  (A/T  ) 

S  cU 

+1.66  log  (  £*)  wts  found  to  agree  quite  well  with  M'  ,  = 

2.4  +  lo*  (A/T4-12  f  6  log.  ( /^  ) .  Moreover,  the  1.6b 

value  was  verified  by  measuring  the  slope  of  (A/r)r-.  versus 
A  (log-log)  for  eit,  .teen  different  events.  No  attempt  was 
made  to  measure  the  same  period  wavelet  at  different  stations. 
The  resulting  spread  in  periods  ranged  from  ;-12  seconds  out 
remained  relatively  _ or.stant  at  a  given  station  for  all 
events.  Next,  average  surface  wave  magnitudes  were  c  „  ..- 
lated  using  M'a  from  aval ’.at  le  :i  mental  .WSSN  stations 
(Table  l).  A  master  ... '.e  »v.u  u  iag  known  yields  ver¬ 
sus  M'  (avg)  (Figux  . '.'..a  f,-  ratios  versus 

known  yi«  Id  were  plotted  for  eight  stations  (Figures  j.2  - 
3.10).  From  the  resulting  nine  master  xurvea,  eleven  unknown 
yields  were  predicted  and  tabulated  in  Table  g 

lone  1 uslons 

A  revised  surface  wave  formula  appears  to  be  borne  out 

by  experimental  data  when  resu’.ts  are  based  on  maximum  A/T 

v 

ratios.  Using  these  ratios,  the  yields  of  eleven  Nevada 
Test  Site  nuclear  events  were  predicted.  The  results  in 
Table  2  indicate  that  the  events  bOXf/A.  RILE  DRIVER,  DURYE/ , 
CHARTREUSE,  and  SCOTCH  are  predicted  wltniu  the  error  bounds 


of  on"  standard  deviation.  However,  HALFbEAK  is  overestimated 
while  GRKKLE i  is  underestimated.  Clearly,  the  majority  of 
Known  yield  event.,  are  predictable;  nothing  can  be  inferred 
fr’oin  the  predictions  of  unknown  yields  until  further  declas¬ 
sification  is  forthcoming.  The  set  of  master  curves  (log-log) 
are  somewhat  inconclusive  below  100  kiloton  where  tne  straight 
line  appears  to  curve  concave  downward.  A  possible  inter¬ 
pretation  would  be  an  Increase  of  seismic  efficiency  with 
decreasing  yield.  That  is,  less  energy  would  be  utilized 
In  inelastic  processes  such  as  rock  vaporization,  plastic 
flow,  and  crushing  and  cracking  rock  material  in  the  near 
source  region.  The  semi-log  master  curve  tends  to  reinforce 
this  hypothesis  by  showing  increased  curvature  Pt  high  yields 
or  less  efficient  seismic  coupling. 

In  the  future,  a  Knowledge  of  more  yields  In  the  40-1000 
kiloton  range  would  greatly  enhance  the  resolution  of  the 
master  curves  especially  near  the  change  of  slope  which  ap¬ 
pears  between  b0-100  Kilotons. 
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Table  2.  Predicted  Yields. 
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4.  Amplitude  Equalization  of  P-wave  Spectra  from 
Underground  Nuclear  Explosions  Recorded  at 
Teleseismic  Distances. 

by  Theron  J.  Bennett 
Introduction 

On  the  basis  of  body  wave  magnitude  the  twelve  largest 
underground  nuclear  explosions  to  date  were  selected.  It 
was  later  found  that  of  these  twelve  only  eight  gave  usable 
records  on  long-period  instruments  at  epicentral  distances 
greater  than  20°.  These  events  along  with  available  informa¬ 
tion  on  times,  coordinates,  and  assigned  body  wave  magnitudes 
are  presented  in  Table  1. 

For  each  event,  the  long-period,  vertical  seismograms 

from  all  the  stations  in  the  U.S.C.Q.S.  World  Wide  Network 

were  scanned  end  those  with  good  signal-to-noise  ratios 

were  selected.  The  restriction  of  epicentral  distances  to 

ranges  greater  than  20°,  in  an  effort  to  reduce  crustal 

complications,  does  not  leave  many  stations  with  good  signal- 

to-noise  ratios  even  for  the  very  large  events  considered 

here.  Those  stations  giving  usable  records  are  listed  for 

each  event  in  Table  1.  These  stations  range  in  distance  up 
o 

to  80  for  some  events  with  a  number  of  South  American  sta¬ 
tions  recording  especially  good  signals  from  the  Nevada  Test 
Site  explosions. 

Method  of  Amplitude  Equalization 
The  spectrum  of  the  recorded  P-wave  signal  can  be  rep¬ 
resented,  following  the  method  of  Ben-Menahem  et  al  (1965), 
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as  the  product  of  the  source  spectrum  and  a  aeries  of  trans¬ 
fer  functions  corresponding  to  various  segments  of  the  trans¬ 
mission  path. 

A  (uj)  *  MJcf^  S  <->  (1) 


where  Hjn,  Hrc>  hM»  and  HSC  rePresent  the  transfer  functions 
of  the  recording  instrument,  receiver  crust,  mantle,  and 
source  crust  respectively. 

The  instrument  response  function  following  Kisslinger 
(1967)  is  given  by: 


(2) 


*  «*  ((,  -«•>*)  ] 


whereof  and  <*>ni  are  the  natural  (angular)  frequencies  of 
the  seismometer  and  galvanometer  respectively,  m  is  a  con¬ 
stant  factor  dependent  on  the  peak  magnification  of  the  system, 
{  and  are  damping  factors,  and  c*  is  the  coupling  factor. 

The  values  for  these  constants  at  stations  in  the  World  Wide 
System  are  provided  by  the  U.S.C.G.S. 

The  transfer  function  for  the  receiver  crust  was  deter¬ 
mined  using  a  computer  program  based  on  the  Haskell-Thompson 
method  for  the  crustal  structures  at  each  station  provided  by 
Teng  (1968)  and  Steinhart  and  Meyer  (1961).  The  response  of 
the  source  crust  for  models  of  the  Nevada  and  Aleutian  test 
sites  were  found  by  applying  the  procedure  developed  by  Fuchs 
(1965)  for  an  explosive  point  source  in  a  layered  crust.  A 
typical  transfer  function  for  the  NTS  crust  model  given  in 


< 
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Table  2  is  presented  In  Figure  4.1  for  a  takeoff  angle  of  35° • 
The  mantle  response  (cf.  Teng  (1968))  is  represented  as 
the  product  of  the  frequency  Independent  geometrical  spread¬ 
ing  and  the  anelastic  attenuation  which  is  a  function  of 
frequency: 


where  ds  is  an  element  of  length  along  the  ray  path,  Q(r)  is 
the  intrinsic  quality  factor  as  a  function  of  depth,  and 
V(r)  is  the  wave  velocity  as  a  function  of  depth. 

The  geometrical  spreading  was  determined  from  an  expres¬ 
sion  given  by  Bullen  ( 1963 ) : 


for  a  surface  focus  source  where  r  is  the  radius  of  the 

o 

earth,  A  is  the  epicentral  distance,  T  the  travel  time, 

and  e0  ls  the  anSle  of  emergence  of  the  ray  at 

the  surface  which  is  given  by  •  c©i“;  RM m. 

The  travel  time  derivatives  were  obtained  by  fitting  a  poly- 
o 

nomial  over  5  segments  of  the  Herrin  P  tables  and  performing 
the  appropriate  differentiations.  This  procedure  is  dis¬ 
cussed  in  greater  detail  in  the  appendix. 


The  integral  in  the  exponent  of  the  anelastic  attenua- 

f  ^  *  - 

tion  factor  of  eq.  3  -  viz.  J  © rr>V(v>  -  was  evaluated 

fSL 

by  numerical  integration  for  the  Q  model  given  in  Table  3  and 


the  Herrin  velocity  model.  (See  appendix  for  details). 
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Signal  Processing 

The  P-wave  signal  for  the  long-period  vertical  record 
was  digitized  at  0.2  sec  intervals  over  a  time  window  of 
approximately  20  sec.  It  is  believed  that  most  of  the  P- 
wave  signal  energy  has  been  included  and  any  other  phases 
such  as  PP  which  may  arrive  within  tne  window  are  small  in 
comparison  to  the  P  signal  strength. 

The  mean  and  linear  trend  were  removed  from  the  data 
and  a  Hamming  window  applied.  The  results  were  then  put 
into  a  Fast  Fourier  Transform  computer  program  to  obtain 
the  amplitude  spectrum,  which  is  defined  according  to  the 
norm: 

FM-  f~  Jtt  (5) 

and  its  transform  pair  as: 

f(t)  •  (6) 

Of  course,  in  the  case  of  discrete  data  these  quantities  take 
on  analogous  forms  with  finite  sums  replacing  the  Integrals. 

Next,  the  effects  of  transmission  path  were  removed  by 
inverting  equation  2  and  the  amplitude  spectra  at  the  source 
determined  for  each  event  at  each  station. 

Observations  and  Conclusions 

A  total  of  75  P-wave  signals  from  the  eight  events  were 
processed  in  this  manner.  A  typical  amplitude  spectrum  for 
each,  corrected  only  for  instrument  response,  is  presented  in 
Figure  4.2.  An  example  of  the  steps  in  the  amplitude- 
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equalization  procedure  is  shown  in  Figure  4.3*  and  the  source 
spectra  for  the  eight  events  of  Figure  4.2  are  '■resented  in 
Figure  4.4.  Amplitudes  in  all  these  figures  are  given  in 
-sec.  The  spectra  show  the  following  features: 

1.  Most  of  the  uncorrected  spectra  exhibit  a  broad, 
overall  peak  between  frequencies  of  0.1  Hz  and  0.9  Hz.  Most 
of  the  other  fluctuations  in  the  spectra  appear  to  correlate 
better  between  events  at  a  single  station  than  between  sta¬ 
tions  for  a  single  event  (cf.  Figure  4.2).  This  would  suggest 
that  these  variations  are  either  due  to  differences  In  crust 
or  upper  mantle  structures  in  the  vicinity  of  the  receivers 

or  possibly  the  result  of  a  radiation  pattern  in  the  source 
region.  A  similar  phenomenon  was  noted  by  Filson  (1970) 
using  shorter  period  P-waves  from  Soviet  explosions  recorded 
at  array  stations. 

2.  Removal  of  the  receiver  crust  seems  to  have  little 
effect  on  the  general  appearance  of  the  spectra  though  the 
overall  level  is  reduced  (cf.  Figure  4.3). 

3.  Correction  for  the  source  crust  has  a  more  drastic 
effect  on  the  shape  of  the  curves.  The  trough  at  the  low 
frequency  end  of  the  spectrum  is  greatly  reduced  and  in 
many  cases  almost  disappears.  Also,  the  level  of  the  curves 
at  frequencies  above  this  notch  is  cut  by  a  factor  of  nearly 
ten  (cf.  Figure  4.3). 

4.  The  most  pronounced  effect  on  both  the  shape  and 
level  of  the  spectra  is  produced  by  the  mantle  correction. 

The  geometrical  spreading  adjustment  changes  the  level  of  the 
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spectrum  at  a  given  station  and  correction  for  anelastlclty 
Introduces  a  linear  trend  in  the  semi-log  plots  of  the  ampli¬ 
tude  spectra  (cf.  Figure  4.3).  More  dramatic  is  the  general 
smoothing  of  the  minor  variations  discussed  under  point  1 
over  the  entire  frequency  band. 

5.  As  one  might  hope,  the  amplitude  spectra  for  a 
given  explosion  look  more  alike  from  station  to  station  after 
application  of  the  equalization  procedure  (cf.  Figure  4.4). 

The  similarity  in  the  general  level  of  the  curves  in  the  band 
0.1  Hz  to  0.8  Hz  seems  noteworthy. 

6.  In  comparing  the  source  spectra  for  different  events 
recorded  at  different  stations  as  in  Figure  4.4,  there  is  no 
obvious  relation  between  the  yield  of  the  explosion  and  the 
level  of  the  spectra.  In  Figure  4.4  it  can  be  observed  that 
throughout  most  of  the  frequency  band  the  amplitudes  for  the 
Greeley  source  spectrum  determined  at  Caracas  are  higher  than 
those  for  the  other  events  depicted  despite  the  fact  that 
the  Greeley  yield  is  only  825  KT  compared  with  yields  greater 
than  1000  KT  for  some  of  the  other  explosions. 

7.  Even  when  the  source  spectra  for  several  explosions 
are  compared  at  a  single  recording  site,  as  in  Figure  4.5*  the 
relation  of  level  and  yield  is  still  not  apparent.  In  this 
case,  the  Halfbeak  explosion  with  a  yield  of  300  KT  shows  up 
with  nearly  the  same  spectral  level  as  the  Boxcar  event  which 
had  a  yield  four  times  as  great. 

8.  In  the  spectra  of  Figure  4.5  minor  peaks  appear  at  the 
love?  end  of  the  frequency  band  for  all  of  the  NTS  explosions 


65 


Neglecting  for  the  moment  Faultless  which  has  some  odd  fea¬ 
tures,  the  spectra  for  the  larger  events  such  as  Benham, 
Boxcar,  and  Greeley  seem  to  exhibit  multiple  peaks  between 
0.1  Hz  and  0.4  Hz  while  the  spectra  of  the  smaller  events 
Halfbeak  and  Commodore  show  a  better  developed,  single  peak 
centered  at  about  0.2  Hz.  Also  the  spectra  for  the  larger 
events  have  minor  troughs  at  a  frequency  of  0.7  Hz  which 
do  not  show  up  for  the  smaller  events.  Similar  features  can 

be  distinguished  on  many  of  the  spectra  at  other  stations. 

o 

9.  At  epicentral  distances  greater  than  about  25  the 
most  notable  features  of  the  spectra  are  the  trough  near 
0.1  Hz  and  the  minor  trough  between  0.7  Hz  and  0.8  Hz. 

It  is  not  known  at  this  time  whether  these  features  or 
any  other  may  be  significant  in  analyzing  near  source  charac¬ 
teristics  of  underground  nuclear  tests. 

APPENDIX.  Mantle  Transfer  Function 

Geometrical  Spreading.  The  most  significant  cause  of 
amplitude  diminution  in  body  waves  is  geometrical  spreading. 
The  theory  of  geometrical  spreading  of  spherical  wave  fronts 
is  concisely  presented  in  Bullen  (1963).  In  a  recent  study 
by  Duda  (1970) calculations  of  geometrical  spreading  based  on 
the  Herrin  P  times  and  Jeffreys -Bullen  S  times  were  presented 
in  conjunction  with  computations  of  body  wave  magnitude  cor¬ 
rection  factors.  In  the  present  research  some  adjustments 
were  made  on  the  procedure  used  by  Duda,  and  these  changes 
along  with  other  results  which  came  up  in  connection  with 
these  refinements  will  be  presented  here. 
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The  geometrical  spreading  can  be  calculated  from  the 
relation: 


G<  • 


_ 

r4%  7,  s  in  A 


1 


-  »?.*  s.'n*e. 


(Al) 


< 


where  r0  is  the  radius  of  the  earth,  A  the  epicentral  dis¬ 
tance,  eQ  the  angle  of  emergence  of  the  ray  at  the  earth's 
surface,  7,  «  }  *7,  •  r«  /v;  >  vo  the  wave  velocity 

at  the  surface,  the  wave  velocity  at  the  source  depth, 
and  T]_  the  radial  distance  from  the  source  to  the  center  of 
the  earth.  Evaluation  of  G(&)  by  eq.  Al  requires  computa¬ 
tion  of  the  second  derivative  of  the  travel  time.  In  the 
calculations  carried  out  by  Duda  the  first  derivatives  of 
the  travel  times  (dT/dA)  were  taken  to  be  the  time  differ¬ 
ences  between  adjacent  epicentral  distances  divided  by  the 
distance  increment,  and  the  second  derivatives  (d2T/d^2) 
were  found  similarly  by  taking  the  differences  between  values 
of  the  first  derivative  at  adjacent  points.  This  corresponds 
to  taking  the  slope  of  a  straight  line  between  the  two  points 
to  be  the  derivative.  While  this  technique,  is  good  as  a 
first  approximation,  it  did  result  in  some  instability  in  the 
spreading  curves  computed  by  Duda,  especially  at  short  dis¬ 
tances  . 

In  the  present  analysis  it  was  decided  that  a  better 
approach  might  be  to  fit  a  polynomial  over  certain  segments 
of  the  travel  times  and  compute  first  and  second  derivatives 
from  that.  Thus  it  was  assumed  that  the  travel  time  could 
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be  represented  as: 

T  *  a.  ♦  a,  &  ♦  + 

and  the  derivatives  as: 

3a  *  •.  ♦  2  «»  L  *  ,SA,+--  +  nft- 

J*T  „r.-a 

«  a  A 

over  some  range  of  epicentral  distances .  The  interval  upon 
which  the  least  squares  fit  was  performed  was  arbitrarily 
chosen  to  be  ten  points  which  corresponds  to  a  5°  segment  of 
the  Herrin  tables.  At  first  a  second  degree  polynomial  was 
fitted  to  the  data,  but  it  was  later  found  that  a  third 
degree  polynomial  fit  gave  more  consistent  results  when  the 
travel  times  were  calculated  from  the  velocity  model. 

In  Figure  4.6  the  computations,  based  on  this  procedure, 
for  (dT/dA  vs  A  )  and(d2T/d62  vs  A  )  are  presented  for  P 
waves  from  a  surface  focus  event.  The  angle  of  incidence  at 
the  surface  as  a  function  of  epicentral  distance,  obtained 
from  the  relation: 


.  -  fjx  ol  T  \ 
lW  \  ro  / 


(A2) 


(where  vQ  was  taken  as  6  km/sec,  r0  as  6371  km,  and  dT/d  A  has 
units  of  sec/radlan),  is  also  shown. 

The  geometrical  spreading  was  then  obtained  using  equa¬ 
tion  A1  and  is  presented  in  Figure  4.7  for  a  surface  focus 
event.  The  amplitude  at  an  epicentral  distance  of  80°  was 


taken  as  a  reference  and  tne  spreading  normalized  on  that 
value  so  that  the  results  could  be  compared  directly  with 
Duda's  curves.  (Duda  had  chosen  this  distance  as  a  refer¬ 
ence  because  it  corresponded  to  a  range  in  which  his  calcu¬ 
lations  were  most  stable.) 

In  comparing  the  spreading  computed  here  with  that  ob¬ 
tained  by  Duda,  the  most  obvious  difference  is  the  improve¬ 
ment  in  stability  as  the  new  curve  is  relatively  smooth  over 
the  entire  epicentral  distance  range.  The  overall  shape 
and  level  of  Duda's  curves  are  maintained. 

Anelastlc  Attenuation.  It  has  been  clear  from  the  earli¬ 
est  days  of  instrumental  seismology  that  the  frequency  con¬ 
cent  of  events  recorded  at  near  distances  differs  signifi¬ 
cantly  from  that  of  the  same  events  recorded  far  away.  It 
is  now  recognized  that  the  mantle  acts  as  a  low-pass  filter 
for  teleseismic  waves  and  this  is  caused  primarily  by  an- 
elastic  behavior  in  the  mantle. 

In  this  research  the  procedure  of  Teng  (1968)  was  used 
to  evaluate  the  effect  of  anelastlc  behavior  on  the  P-wave 
spectra.  Teng  had  used  the  CIT11A  velocity  model  in  his 
analysis  while  in  this  study  the  Herrin  model  was  used;  it 
was  therefore  necessary  to  make  some  adjustments  in  the 
Q(r)  model  adopted  by  Teng. 

In  equation  3  it  was  noted  that  the  effect  of  anelastlc 
attenuation  could  be  represented  Dy  a  factor  of  the  form: 
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The  expression: 

t*  *  *■  / G/rt-M  (M) 

pk 

where  the  *<  refers  to  P-waves,  was  evaluated  numerically 

for  several  of  the  Q  models  suggested  by  Teng  and  a  quantity 

called  the  differential  attenuation  as  a  function  of  A, 

o 

normalized  to  an  epicentral  distance  of  50  , 

SAC*)*  t*Ca&-  t *(«••>  (A5) 

was  calculated  for  ray  paths  from  an  event  with  focal  depth 
of  550  km.  This  quantity  was  then  compared  to  observed 
differential  attenuations  given  by  Teng  for  two  deep  focus 
earthquakes  in  South  America.  As  might  be  expected,  since 
the  velocities  of  the  CIT  11A  model  and  the  Herrin  model  do 
not  differ  drastically,  the  adjustment  in  the  Q(r)  structure 
required  to  obtain  a  fit  to  the  same  attenuation  data  was 
not  large.  The  best  fit  (cf.  Figure  4.8)  seemed  to  corre¬ 
spond  to  the  model  which  Teng  had  suggested  as  an  alternate . 
This  model  is  shown  in  Table  3. 

Using  this  Q  model  t*  as  a  function  of  epicentral  dis¬ 
tance  was  calculated  for  a  surface  focus  event.  The  results 
are  shown  in  Figure  4.9.  Then  the  expression 

e*r  {•  f  t*} 

giving  the  anelastic  part  of  the  mantle  transfer  function 
can  be  directly  evaluated.  To  show  how  this  affects  the 
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spectra  of  P-waves  at  different  epicentral  distances,  the 
mantle  transfer  function  was  computed  taking  into  account 
only  the  anelastic  behavior  for  distances  of  5°,  10°,  20°, 
and  50°.  The  results  are  presented  in  Figure  4.10.  The 
low  pass  filter  character  is  evident. 
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Table  1 


Event 

Date 

Time 

(GMT)  mh 

Location 

Stations 

Analyzed 

Benham 

(Ben) 

12/19/68 

16:30:00.0  6.3 

37.232N 

116.474W 

ARE, BOG, CAR, 
FLO, NNA, OXF, 
QUI,SHA 

Boxcar 

(Box) 

4/26/68 

15:00:00.0  6.3 

37.295N 

116.456W 

AAM, ARE, BOG, 
CAR, COL,  FLO, 
OXF,  QUI,  SHA 

Commodore 

(Com) 

5/20/67 

15:  00:  00.0 

37.130N 

II6.O63W 

FLO,LPS,  OXF, 
SHA 

Faultless 

(FLT) 

1/19/68 

18: 15:  00.0 

38.625N 

116.219W 

AAM, ARE, ATL, 
BLA,BOG, CAR, 
COL,  FLO, NNA, 
OXF,  QUI, SHA, 
UNM, KRK 

Greeley 

(ORE) 

12/20/66 

15:30:00.1  6.05-6.29 

37.302N 

116.408W 

AAM,  ARE,  ATL, 
BLA,BOG,CAR, 
COL, FLO,LPS, 
NNA, OXF, QUI, 
SHA 

Ha If beak 
(HLF) 

6/30/66 

22:15:00.0  6.02-6.04 

37.316N 
116. 299 Vv 

AAM,BLA,LPS, 
OXF,  SHA 

Jorum 

(JOR) 

9/16/69 

14:30:00.0  6.2 

37.314N 

116.460W 

ATL,BLA,COL, 

FLO,SCP,SJG, 

WES 

Milrow 

(MIL) 

10/2/69 

22:06:00.0  6.5 

51.417N 

179.182E 

ALQ,BLA,COL, 
DUO,FLO,HNR, 
JCT,LON, NDI, 
OXF, QUE,SHA, 
SHL,TUC 
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Table  2 


NTS  Crust 

Model 

Layer  Thickness 
D(km) 

P-Vel 
«*(kro/sec ) 

. .  Kmskehi 

THmBSi 

2.14 

3.24 

1.86 

2.20 

0.36 

4.40 

2.50 

2.40 

2.50 

5.10 

2.60 

2.50 

20.00 

6.10 

3.50 

2.70 

10.00 

7.00 

4.00 

3.00 

• 

8.00 

4.60 

3.30 

Table  3 


Mantle  Q 

Model 

Depth  (km) 

Q 

Depth  (km) 

Q 

Depth  (km) 

Q 

0 

450 

1500 

1000 

287* 

100 

36 

60 

1700 

1000 

50 

75 

1875 

750 

70 

75 

2000 

500 

150 

75 

2200 

500 

530 

75 

2400 

500 

750 

100 

2500 

300 

950 

250 

3600 

250 

1100 

450 

2700 

200 

1300 

550 

2800 

120 

RDINS 


7^ 


FACTOR 


(DE6) 
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Transmission  In  Crustal  Layers 
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5.  The  Effect  of  a  Dipping  Layer  on  P  Wave  Transmission. 

Albert  Rogers 

Crustal  Inversion  utilizing  the  amplitude  spectra  of 
P  or  S  body  waves  makes  possible  the  determination  of  the 
crustal  parameters,  depth  and  velocity  contrast,  beneath  a 
single  station  with  data  from  that  station  alone.  It  is, 
therefore,  not  only  a  very  convenient  way  to  measure  qualities 
that  are  needed  In  other  seismologies 1  studies,  but  it  also 
provides  a  means  of  adding  to  basic  knowledge  about  the  crust. 
The  method  la  particularly  useful  in  remote  areas  of  the 
earth  where  refraction  and  reflection  surveys  are  difficult 
to  conduct. 

A  number  of  researchers,  Phlnney  (l),  Fernandez  (2,3*4) 
Leblanc  (5)*  Hasegawa  (6,7)*  Ellis  and  Bcsham  (6),  Castano 
(9),  Ibrahim  (10),  and  Utau  (11),  have  attempted  to  correlate 
the  theoretical  results  predicted  by  the  Haskell -Thoms on  formu 
lstlon  applied  to  body  waves  with  body  wave  observed  spectra. 
In  some  cases  the  lack  of  a  correlation  may  Imply  the  exis¬ 
tence  of  a  more  complex  model.  Such  complexities  as  aniso¬ 
tropy,  snelsstlclty,  and  heterogeneity  within  the  crustal 
system  may  significantly  affect  the  amplitude  response  of 
body  waves,  Introducing  errors  into  the  inversion  process. 

The  distortion  of  the  transfer  function  may  be  great  enough 
to  make  Inversion  Impossible.  The  introduction  of  dip  into 
one  or  more  of  the  layer*  Is  one  such  complication. 


In  this  study  we  examine  the  changes  produced  in  the  4( 

transfer  function  for  P  waves  passing  through  a  dipping 
layer.  Two  aspects  are  considered:  the  errors  due  to  dip 
introduced  into  the  crustal  depth  determination  when  the 
layer  is  erroneously  assumed  to  be  horizontal;  and,  the 
applicability  of  a  ray  theory  embodying  dip  to  improve  the 
accuracy  of  the  depth  determination  as  well  as  to  find  the 
dip  beneath  the  station. 

We  assume  a  single  dipping  layer,  that  is,  a  wedge, 
overlying  a  half-space.  Plane  P  waves  are  incident  from 
below  the  layer,  and  the  diffracted  waves  that  are  generated 
at  the  vertex  of  the  wedge  are  ignored.  The  difficulties 
associated  with  finding  the  normal  mode  solution  make  a 
first  approach  involving  analogue  seismic  models  and  ray 
theory  very  attractive. 

The  ray  theory  we  have  developed  is  a  computational 
scheme  for  summing  the  plane  wave  P  and  SV  potentials  subject 
to  the  boundary  conditions  expressed  by  Knott's  equations 
and  Snell's  law.  We  have  used  the  theory  to  determine  the 
vertical  and  horizontal  displacements  for  the  multiple  re¬ 
flection  solution  for  a  dipping  layer.  In  the  following 
discussions  the  model  study  transfer  function  is  defined  as 
the  ratio  of  the  vertical  displacement  at  the  top  of  the  layer 
to  the  total  incident  P  displacement  at  the  base  of  the  layer. 

For  the  crustal  study  the  transfer  function  is  defined  as 
the  ratio  of  the  vertical  transfer  function  to  the  hori¬ 
zontal  transfer  function. 
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Matching  of  the  observational  transfer  function  to  the 
theoretical  may  be  accomplished  by  eye  or  analytically. 
When  matching  by  eye  the  theoretical  transfer  function  is 
plotted  on  log-log  paper  in  dimensionless  frequency  Hf/v, 
where  H  is  the  depth  beneath  the  station,  v  is  the  layer 
velocity,  and  f  is  frequency.  One  such  curve  represents 
the  theoretical  transfer  function  for  any  layer  depth  and 
velocity,  for  a  given  velocity  contrast,  angle  of  incidence, 
and  angle  of  dip.  The  observational  curve  is  plotted  on 
log-log  paper  of  the  same  cycle  size.  The  curves  are  then 
superimposed  and  translated  without  rotation  until  a  fit  is 
observed.  The  depth  is  determined  by  reading  from  the  axes 
any  two  corresponding  values  of  Hf/v  and  f,  and  assuming  a 
value  for  the  layer  velocity. 

A  more  objective  method  was  developed  using  the  corre¬ 
lation  coefficient  and  the  mean  square  deviation  defined  by 
(1): 

Correlation  coeff  a  cc.  » 

(i) 


Mean  sq.  dev.  s  ^  2  -jp  ^ 

OtSI  f 


where  y^  equals  the  theoretical  transfer  function  and  yj  equals 
the  observational  transfer  function.  The  correlation  coeffi¬ 
cient  is  computed  about  the  mean  y  ,  and  normalized  by  the 
standard  deviation  so  that  it  is  primarily  a  measure  of  the 
degree  of  equality  in  shape  of  the  two  variables  y^^  and  yj. 

The  mean  square  deviation  chiefly  determines  the  degree  of 


^quality  In  .imp  i  1  tud<-  :<iwecn  me  two  urves. 

Tne  r:iea:.  deviation  was  computed  after  Tilting 

the  two  transfer  fjn.-tlons  uy  eye  and  norrna llzlng  them  at 
orreapondin^  values  In  tne  range  of  tne  peak  seismograph 
.njgnlf  1  -a  t  lo:i . 

In  order  to  analytically  snlft  the  two  transfer  func¬ 
tions  past  one  another  in  the  search  for  the  oest  fit,  they 
were  linearly  Interpolated  at  intervals  of  .01  log  frequency 
arid  log  dimensionless  frequency.  The  corrections  were  oe- 
termined  for  earthquake  data  for  lag  values  measured  In 
kilome  ters  between  13  km  and  90  km.  This  Is  less  than  1  ‘yjL 
of  the  total  lag  measured  from  the  zero  lag  position.  The 
model  data  lag  values  varied  between  16  and  30  mm,  which  Is 
less  than  +  lt>  of  the  total  lag.  The  number  of  Independent 
points  used  In  the  region  of  maximum  fit  varied  from  40  to  30 
for  the  earthquake  data  and  was  equal  to  23  for  the  model 
djta . 

In  the  evaluation  of  the  ray  theory  we  have  studied  the 
effects  on  the  earthquake  transfer  function  of  variations  of 
velocity  contrast,  depth,  angle  of  emergence,  angle  of  dip, 
and  spectral  resolution.  The  transfer  functions  have  been 
evaluated  for  a  Poisson's  ratio  of  .23  and  a  density  ratio 
of  .63. 

Figure  (3.1a)  shows  the  earthquake  transfer  function 
for  velocity  ratios  of  .9,  .6,  and  .7.  It  Is  apparent  that 
Increasing  the  contrast  Increases  the  level  of  the  peaks 
but  does  not  appreciably  alter  the  shape  of  the  curves.  What 


bb 

differences  are  present  are  due  to  the  fact  that  the  curves 
have  unavoidably  been  s soothed  in  a  slightly  different  manner. 
The  nature  of  the  prograa  for  obtaining  the  curves  also  makes 
it  inconvenient  to  obtain  values  at  exactly  the  same  intervals 
of  dlaensionlesB  frequency.  These  transfer  functions  were 
evaluated  for  10°  of  dip  and  60°  angle  of  emergence. 

Figure  (5.1b )  shows  the  transfer  function  for  two  values 

o  o 

of  depth  with  a  dip  10  and  an  angle  of  emergence  126  •  It 

is  clear  that  changes  in  depth  produce  no  changes  in  the  shape 

or  level  of  the  transfer  function.  Any  changes  which  are 

present  can  be  explained  by  the  argument  given  above.  Likewise, 

we  could  expect  there  to  be  no  change  in  the  transfer  function 

level  or  shape  had  only  the  layer  velocity  been  changed  with 

contrast  velocity. 

Varlationa  of  the  angle  of  Incidence  reveal  that  the 
level  of  the  transfer  function  lncreaaea  without  major  charac¬ 
ter  change  aa  the  angle  of  incidence  varies  from  50°  to  fcO°. 
There  la  a  major  character  change  between  b0°  and  100°  ae  wave 
propagation  goes  from  up  to  down  dip.  As  the  angle  of  inci¬ 
dence  varies  from  100°  to  130°  the  level  of  the  transfer 
function  decreases  without  major  character  change. 

We  used  the  correlation  coefficient  to  determine  the 
theoretical  error  resulting  from  fitting  a  0°  dip  theoretical 
transfer  function  to  0°,  5°»  10°,  15°,  20°,  2^°.  Equivalently , 

transfer  functions  plotted  in  frequency  were  shifted  peat  th« 

«o 

0  transfer  function  plotted  in  dimensionless  frequency  and 
tha  correlation  and  mean  square  deviation  were  determined. 


T.'i*-  .  >  1  .*<- .  »:•  1  v.  r.  I  r.  ?:itlc  .  for  the  up  a.-yd  down  dip  :;•»  e» 

u 

For  »<v<m  If.  l  !«•*.•  a*.  <  *.  }»c  orrvlatlors  ©efficient  r  oe  & 

'  fu'uj/f.  i  -rs  1  rilr. .  v.  u  There  *<*a  aes  to  t<«-  any  svaourat  le  fit 

•  a  m.  lip  v.i r l ••  fro*  rv  °  to  *Jio  fit”  lndicalcu  that 

iii  v.iia«a  of  tfw  urrelatlon  flu  turtle  ii round  zero  or  eaw 

iouvr  venue  wttn  no  well  defined  aa  <  l  must .  The  neon  a  qua  re 

o  o 

deviation  Increases  rapidly  In  the  region  0  to  10  .  la  con- 
o  ^ 

slant  near  it  arud  i  .  .  and  be-:  1  :.»>  to  Increase  s^atn  beyond 
1  The  depth  value#  arc  determined  at  point#  of  Maxim** 

correlation,  and  it  la  evident  that  error  In  tnc  depth  deter¬ 
mination  la  lean  than  'it  In  all  caeca  and  lean  than  2. *>4  for 

O  o 

dip#  lean  than  or  equal  to  1C  .  Kor  waves  ln<  Ident  at  1^0 

tne  correlation  coefficient  tfoeo  tnrourfh  a  alnlwun  or  region 
o  o 

of  no  fit  at  1  .  and  a^al.  at  ?'>  .  The  mean  square  deviation 

Increaaco  continuous  1/  an  the  dip  Increases.  The  error  In  the 

value  of  the  depth  varies  fro*  ?i  to  74.  This  Indicates  that 

In  mos t  Instances  the  depth  of  the  layer  *a y  he  date rained 

within  74  by  Matching  to  the  transfer  function  for  a  horizon- 

o 

tal  layer  even  though  the  la>er  May  have  dips  up  to  K  . 
However,  the  decree  of  correlation  May  be  snail  or  zero  in 
some  cases.  It  is  this  fact  that  nay  permit  the  determina¬ 
tion  of  the  dip  as  well  as  the  depth. 

It  is  useful  to  study  the  value  of  the  correlation  co¬ 
efficient  fron  a  statistical  viewpoint.  Vc  nay  test  whether 
the  value  of  rj^  differs  significant ly  fro*  zero  or  whether 
two  rjj  value 8  differ  slr;nif Icantly  fro*  each  other.  *e  may 
also  determine  confidence  limits  of  the  correlation  in  the 


VO 


population,  If  w  tett  tha  null  hypothesis  that  the 
theoretical  tranafar  function*  compared  above  are  uncorrelated, 
that  la  drawn  from  a  normal  bivariate  population  with  ^r*  -  0. 

against  the  alternative  hypothesis  M  that  t hey  are  positively 
correlated,  It  la  found  that  a  positive  correlation  exists  In 
all  eases  except  for  the  instances  of  no  fit.  If  the  value 
of  >  .214  (this  number  applies  only  to  these  theoretical 

curves),  H  la  accepted  at  tha  level  of  slgnlf lcance .  * 

value  of  rjj  this  large  due  to  chance  would  occur  only  once 
In  twenty  times  if  the  sample  la  drawn  frow  a  population  with 


It  la  wore  Informative,  perhaps,  to  tv  t  If  tha  measured 
values  of  the  correlation  are  taken  frow  a  population  of  high 
correlation,  aaj  ^r^  -  .9*  The  9Vd  confidence  limits  given 
In  Table  1  Indicate  thet  only  tha  *>°  correlations  (Including 
the  0°  cases)  for  both  angles  of  emergence  could  have  been 
taken  from  a  population  of  high  correlation.  In  all  other 
cases  we  would  be  forced  to  eccept  the  alternative  hypothesis 
of  low  correlation.  The  Haskell -Thomson  model  yields  statis¬ 
tically  algnlf leant  correlations  when  dipt  less  than  5°  are 
present.  9ot  dipt  greater  than  *>0  their  model  does  not  give 
significant  correlations,  and  in  soma  eaaea  the  correlation 
la  low  or  aero.  Theoretically,  even  in  cases  of  low  correla¬ 
tion,  the  depth  may  be  determined  within  Til. 

In  the  theoretical  studies  of  the  spectral  resolution 
we  measured  the  effect  on  the  transfer  function  of  applying 
a  Hamming  window  of  60,  *0,  and  20  sec  ler*th  to  theoretical 


sr l saoff.raas  .  airylowa  as  short  as  i*c  sec  have  no  effect  oei  the 
accvira-y  of  the  depth  determination  for  •  k*  crust.  How¬ 
ever.  i»  short  window  may  smooth  the  transfer  function  enough 
to  wane  a  dip  deterslfiSt  Ion  Impossible. 

The  analogue  rode  1  study  was  conducted  In  an  effort  to 
Investigate  under  am  rolled  laboratory  conditions  the  see  ora:/ 
of  the  -rusisl  Inversion  process.  It  siso  provided  s  test  of 
the  ray  theory.  The  experiments  were  designed  for  the  measure¬ 
ment  of  the  vert  ice  1  trsnsfer  function  for  three  values  of 
o  o  o 

dip.  n  ,  7 . .  ic  .  and  for  vertically  Incident  eaves.  To 
what  extent  the  model  results  apply  to  crustal  transfer  func¬ 
tions  Is  unknown;  however,  subsequent  analysis  reveals  that 
conclusions  baaed  on  the  modal  results  are  entirely  comparable 
to  the  conclusions  reached  using  the  theoretical  crustal 
transfer  functions. 

For  the  experiments  wa  used  a  two  dimensional  model. 

Figure  (*>.2),  with  a  1/e  Inch  title  H  aluminum  alloy  ha  If- space 
and  a  Plexiglas  layer.  The  model  was  actually  constructed 
with  the  solid-solid  Interface  horizontal  and  the  dip  In  the 
free  edge  of  the  layer.  The  depth  from  the  receiver  to  the 
half-space,  measured  along  the  normal  to  the  free  surface, 
was  constant  for  all  three  cases.  For  the  free  surface  to 
represent  a  horizontal  plane,  at  In  the  reel  earth,  the  re¬ 
ceiver  was  displaced  from  the  model  vertical  by  en  angle 
equal  to  .he  angle  of  dip.  These  two  steps  insured  that  the 
amplitude  response  that  was  observed  reflected  only  changes 
in  the  sngle  of  dip.  The  lower  receiver  measured  the  source 


function  end  mss  the  mm  distance  f pom  the  source  ss  the 
source  •«§  fro*  the  base  of  the  layer  beneath  the  receiver. 

The  two  cent lneter  layer  thickness  beneath  the  receiver 
was  chosen  ss  s  compromise.  Zf  the  layer  depth  were  Much 
shallower,  the  ms Jlor  spectral  peaks  and  troughs  wouH  become 
widely  spaced  snd  .'.hi  ft  to  high  frequencies  (based  on  the 
formula  for  normal  incidence  to  a  horizontal  layer).  The 
source  function  would  llluMlnate  few  or  none  of  the  peaks 
snd  troughs  In  such  a  spectra*.  On  the  ©the  hand.  If  the 
layer  depth  were  Made  Much  greater,  the  peaks  and  troughs 
would  beeches  ao  closely  spaced  that  they  could  not  be  re¬ 
solved  by  the  available  spectral  resolution.  In  addition, 
the  two-way  travel  tlwe  increases  snd  few  Multiples  would  be 
observed  before  interfering  phases  fraw  the  sides  of  the 
Model  arrive. 

Capacitive  receivers  ware  used,  end  sliver  scetyllde 
explosive  charges  served  ss  the  energy  source.  When  properly 
formed,  these  charges  yielded  e  slMple  broad -band  signal, 

5  kllx  to  300  kHt,  that  adequately  1  11um  Ins  ted  the  layer  sMpl  1  • 
tude  spectruM. 

A  comparison  between  the  theoretical  snd  observed  model 

transfer  function  is  plotted  In  figure  (5.3).  The  degree  of 

fit  appears  "best"  In  the  snd  10°  esses.  The  7.^°  results 

are  not  ss  significant  due  to  the  interference  of  antisymmetric 

o  o 

modes.  This  problem  was  ove rooms  for  the  5  *nd  10  experi¬ 
ments.  The  regions  of  significant  deviation  are  below  kHz 
ss  evidenced  by  sn  Increase  in  saplitude  and  shift  to  lower 


frequency  in  relation  to  lh»  lit* orotic*!.  HoMtvtr,  the  ray 
theory  la  not  valid  in  this  frequency  band,  since  tha  mve> 
ler^tna  ere  2  to  l)  times  tha  layar  thickness.  Above  160  kHz 
there  is  a  distinct  frequency  dapandant  shift  of  tha  obaarvad 
peaks  and  troughs  towards  fraquanciaa  lower  than  tha  thaorat V 
cal.  This  la  accounted  for  by  tha  occurrence  of  geometrical 
dlsparalon  and  tna  breakdown  of  two -dimensionality  in  tha 
model.  Tha  condition  for  daatructiva  lntarfaranca  for  normal 
lncldanca  to  a  horizontal  layar  is  divan  by, 

r  vah 

and,  t hare fore,  a  change  In  valoclty  producaa  a  proportional 
change  In  tha  frequancy  of  lntarfaranca . 

4*  <  4»An 

o* 

o 

Comparing  tha  obaarvad  and  experimental  for  5  ravaala  that 
tha  trough  at  230  kHz  (theoretical)  appears  at  211  to  216  kHz, 
a  proportional  decrease  in  percentages  of  6-60.  Tha  phase 
valoclty  of  a  1/8"  thick  Plexlflaa  layar  at  230  kHz  la  deter¬ 
mined  from  tha  aquation  (aftar  Tolstoy  and  Usd  In  (12)  )i 

TmJH'  -<*•*!*  WAJ  [«  *  <-VyO* 


Tha  taro  frequancy  limit  of  tha  phase  valoclty  la 

C •  *  tVa /  <!  a^i»l 


The  principle  root  of  the  transcendents  1  elution  us  In; 
the  aodel  constants  and  230  kHz  lac  -  i.‘,  c  am/ /•  z .  This 
represents  s  decrease  in  tne  velocity  of  tt. 

The  correlation  coefficient  technique  was  employed  to 
analytically  Measure  tne  error  In  the  depth  determine  t Ion 
as  well  as  to  study  the  feasibility  j:  -s-raa  .,r  Ing  the  dip. 
Table  (2)  shows  the  results  of  rossp-arirv  va  t.  of  the  experi¬ 
mental  spectra  to  their  corresponding  theoretical  spectra. 
Only  the  valuea  up  to  160  kHz  were  used  In  order  to  exclude 

the  effect  of  geometrical  dispersion.  As  was  observed  by 

o  o 

eye,  the  degree  of  fit  la  better  In  the  *>  and  10  spectra 
according  to  the  Measured  valuea  of  the  correlation  coeffi¬ 
cient  and  the  mean  square  deviation.  The  Measured  depths 
of  the  layers  beneath  the  receiver  wire  20.0  am.  20. }  mm, 
and  20.5  am.  for  the  t>°,  7.3°,  and  10°  lsyera  respectively. 
The  measurement*  Include  the  epoxy  bonding  layer  aa  well  as 
the  foil  used  as  one  piste  of  the  receiver  pickup.  The 
values  of  depth  determined  using  the  correlation  coefficient 
are  In  error  by  about  ♦  1JK.  An  error  of  this  size  may  be 
made  In  (maturing  tne  depth;  In  addition,  the  correlation 
was  computed  In  Intervals  of  ,**6  mm,  and.  therefore,  the 
depth  aey  be  In  error  by  *  .23  mm  or  approxlaetely  ♦  14. 

It  la  evident  that  the  depth  of  a  layer  may  be  determined 
to  e  high  degree  of  accuracy  In  a  controlled  experiment  even 
though  the  correlation  may  be  eamii.  This  conclusion  Is  In 
sgreeaent  with  previous  theoretical  results. 


f') 


£  act  of  itfk&SM  earthquake  records  were  chosen  and 
arw  iyccd  in  order  to  determine  the  practicability  of  apply- 
in*  these  results  to  the  real  earth.  Long  period  records 
were  chosen  even  though  portions  of  the  spectral  band  violate 
the  regions  of  validity  of  the  ray  theory  (wavelength  such 
longer  than  the  layer  tnicknesa).  The  frequencies  uaed  were 
.012*5  cpa.  to  .o  *po  which  repreaent  waveler^tha  varying  froa 
10  to  1/*)  the  layer  thickness.  Table  (3)  givea  all  the  per¬ 
tinent  earthquake  atatistlca. 

The  theoretical  transfer  functiona  were  evaluated  uaing 
tne  required  angle  of  incidence,  a  fixed  velocity  con treat 
of  .  b,  and  a  fixed  denally  contraat  of  .63.  Table  (4) 
contalna  the  reaulta  of  coopering  theoretical  to  obaervad 
tranafer  functiona  for  dlpa  of  0°  through  30°  in  5°  Incre¬ 
ments  for  two  earthquakes  recorded  at  Antofagaata .  Two 
features  clearly  stand  out.  first,  the  corrolstlon  o ©effi¬ 
cients  have  double  valued  aoxlas,  and  aecorr  ,  the  correlation 
goaa  through  a  minimum  at  10°.  The  double  values  are  due 
to  the  "periodic"  nature  of  the  tranafer  function.  Coopering 
the  two  regions  of  fit  by  eye  indicates  tha  fit  for  tha 
smaller  deptn  to  be  beat.  Geologically,  tha  smaller  depth 
aeeme  wore  likely  since  Antofagaate  is  on  tha  Chilean  coast 
about  120  km  f row  the  trench  axis.  However,  there  la  lost 
evidence  suggesting  tha  isrgtr  depth  baaed  on  ref.'tctlon 
and  gravity  surveys.  This  evidence  will  be  discussed  fur¬ 
ther,  subsequently .  The  fact  that  tha  value  of  the  correla¬ 
tion  coefficient  goea  through  a  minimum  at  10°  also  noted 
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In  Table  (1)  where  two  theoretical  transfer  functions 
were  compared.  Those  functions  are  also  the  theoretical 
curves  for  Antofagasta.  This  Indicates  that  the  increase 
In  the  value  of  the  correlation  above  10°  cannot  be  asso¬ 
ciated  with  actual  dips  that  large.  Furthermore,  the  con¬ 
fidence  Intervals  Indicate  that  for  the  earthquake  of 
o  o 

17/01/67,  the  0  and  9  spectra  represent  a  significantly 
better  fit  at  the  95£  level  than  do  any  other  cases.  The 
confidence  Interval  also  Indicate  that  va  cannot  statistically 
differentiate  between  0°  and  5°  of  dip;  however,  we  may  accept 
the  valuea  aa  "Indicators"  of  the  true  dip.  The  5°  case 
appears  to  be  slightly  more  significant  when  both  the  cor¬ 
relation  and  mean  square  deviation  are  taken  Into  account 
for  both  earthquakes.  This  result  agrees  with  the  "best" 
fit  by  eye.  The  vslues  of  the  depth  that  are  obtained  for 
5°  of  dip  are  47.2  km  and  46.0  km.  The  mean  Is  46.1  with  a 
atandard  deviation  of  ♦  1.5  km.  The  final  fit  that  Is  ob¬ 
tained  for  both  earthquakes  Is  shown  In  Figure  (5.4). 

The  results  for  earthquakes  recorded  at  Nana,  La  Paz, 
and  Arequlpa  are  given  In  Table  (5).  The  correlation  co¬ 
efficients  and  the  95£  confidence  Intervals  for  NNA  reveal 
that  the  0°  and  5°  theoretical  carves  do  not  correlate  posi¬ 
tively  with  the  observational  spectrum.  Positive  correla- 

o  o  o 

tlona  do  exist  for  10  ,  15  »  and  20  of  dip  with  the  maximum 
o 

occurring  at  15  »  although  the  maximum  Is  not  significantly 
different  than  the  other  values  at  the  95$  level.  Taking 
Into  account  both  the  correlation  coefficient  and  the  mean 
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square  deviation,  the  best  fit  is  for  a  layer  dipping  15 
SE  and  a  depth  of  74.7  km.  Figure  (5.5)  shows  the  chosen 
fit. 

Positive  correlations  are  found  for  LPB;  however,  the 
degree  of  visual  correlation  shown  in  Figure  (5.5)  is  poor. 
In  the  case  of  ARE  there  are  no  positive  correlations  and 
the  visual  fit  is  also  poor.  Hence  no  conclusions  regarding 
these  data  will  be  drawn.  Both  of  these  stations  are  at 
high  elevations,  and  it  is  likely  that  the  crustal  struc¬ 
ture  beneath  them  is  complex. 

The  crustal  depth  determined  at  Antofagasta  is  in  good 
agreement  with  results  of  refraction  data  slightly  to  the 
northeast.  Aldrich  et  al  (13)  give  a  value  of  46  ion  based 
on  a  one  layer  crust,  with  velocity  6.0  km/sec.  There 
are  indications  of  a  two  layer  crust  with  a  second  layer 
velocity  of  7.0  km/sec  and  total  crustal  thickness  of  56  km. 
Our  data  possibly  have  two  interpretations  as  well,  since  a 
fit  was  also  obtained  for  a  57-3  km  crust.  Study  of  two 
layer  transfer  functions  (Fernandez  (14))  indicate  that  for 
some  velocity  ratio  and  depth  combinations  it  is  possible 
that  the  transfer  function  became  "doubly  periodic"  in  re¬ 
lation  to  the  one  layer  transfer  function.  That  is,  for 
each  peak  in  the  one  layered  transfer  function,  there  are 
two  in  the  two  layered  case.  Therefore,  in  fitting  a  one 
layered  theoretical  transfer  function  to  a  two  layered  ob¬ 
servational  spectrum,  two  matches  may  actually  be  possible. 
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Gravity  data  interpreted  by  Lomnitz  (15)  near  Anto¬ 
fagasta  also  support  a  larger  value,  near  55-60  km. 

The  dip  of  the  Moho  near  both  the  Peruvian  and  Chilean 
o  o 

coasts  is  between  5  and  10  using  data  and  interpretation 
of  Fisher  and  Raitt  (l6),  Woolard  (17),  and  Tatel  and  Tuve 
(18).  Lomnitz  ( 15 )  north-south  gravity  profile  in  Chile 
indicates  a  crust  dipping  gently  (approx. 1/2°)  to  the  north 
leading  to  crustal  depths  of  70  km  beneath  northern  Chile 
and  southern  Peru.  Our  75  km  depth  determination  at  Nana 
indicates  that  a  deep  crust  extends  into  northern  Peru  as 
well.  Geologically,  we  may  expect  a  deep  crust  at  Nana, 
since  it  is  located  at  the  foot  of  the  Andes  about  240  km 
from  the  trench  axis  -  twice  as  far  as  Antofagasta  is  from 
the  trench  axis. 
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6.  Geological  Factors  Affecting  Triggering  and  a 
Model  of  Aftershock  Production, 
by  Carl  Klssllnger 

It  is  known  that  when  an  explosion  is  fired  in  a  solid 
medium  under  stress,  the  resulting  seismic  signal  is  altered 
from  that  produced  by  a  similar  explosion  in  an  unstressed 
medium.  With  reference  to  the  specific  problem  of  monitor¬ 
ing  underground  explosions,  the  important  question  is  whether 
this  alteration  of  the  signal  can  make  the  detection,  loca¬ 
tion,  or  identification  of  the  event  more  difficult.  The 
principal  concern  is  whether  the  concurrent  release  of  tec¬ 
tonic  strain  energy  at  the  time  of  the  explosion  can  suffi¬ 
ciently  alter  the  seismic  signal  to  make  application  of 
identification  criteria  less  certain.  A  closely  related 
question  is  whether  it  might  be  possible  to  Intentionally 
place  an  explosion  in  a  tectonic  environment  in  which  a  large 
release  of  tectonic  strain  by  faulting  would  occur,  l.e.  a 
large  earthquake  would  be  triggered  intentionally,  as  an 
evasion  technique. 

An  ambient  field  stress  field  will  affect  the  seismic 
signal  through  several  mechanisms.  Of  these,  the  least 
important  for  the  present  problem  is  the  introduction  of  a 
small  amount  of  velocity  anisotropy  because  of  the  stress 
dependence  of  the  elastic  moduli.  More  Important  is  the 
control  exerted  by  the  ambient  stress  on  the  generation  of 
tensile  fractures.  In  a  homogeneous  isotropic  unstressed 


material,  the  tensile  fractures  that  mark  the  outer  limits 
of  non-elastic  behavior  will  be  distributed  more  or  less 
uniformly  around  the  shot  cavity,  with  essentially  the  same 
length,  so  that  the  basic  spherical  symmetry  of  the  source 
is  not  markedly  disturbed.  If,  however,  an  appreciable 
prestress  is  present,  tensile  cracks  in  a  few  directions  will 
grow  to  considerable  length.  Prestress  is  not  the  only  con¬ 
dition  that  can  result  in  preferential  fracturing.  Planes 
of  weakness  resulting  from  depositlonal  processes  will  have 
a  similar  effect,  especially  in  surficial  materials. 

A  propagating  tensile  fracture  is  a  source  of  seismic 
waves  and  these  will  be  superimposed  on  the  primary  signal 
from  the  explosion.  This  is  one  mechanism  by  which  S  waves 
are  generated  directly  by  the  explosion  (Kisslinger,  et  al, 
1961).  The  energy  for  these  waves  is  derived  from  the  expand¬ 
ing  explosion  products  and  the  role  of  the  ambient  stress 
is  primarily  that  of  determining  the  direction  of  crack 
propagation. 

Finally,  prestress  will  result  in  the  production  of  a 
seismic  signal  due  to  release  of  stored  elastic  strain 
energy  in  response  to  the  shot .  One  way  in  which  this  will 
occur  is  as  a  result  of  the  sudden  creation  of  an  enlarged 
cavity  in  the  medium.  Another,  potentially  more  dramatic, 
is  the  occurrence  of  slip  on  fault  surfaces  in  response  to 
the  explosive  loading. 

It  is  this  last  effect,  the  triggering  of  a  tectonic 
earthquake  by  the  explosion,  that  is  most  likely  to  cause 


difficulty  in  applying  discrimination  criteria.  However, 
a  concurrent  earthquake  can  cause  difficulty  only  if  It 
generates  seismic  waves  comparable  in  energy  to  those  pro¬ 
duced  by  the  explosion  itself.  An  earthquake  with  magnitude 
equal  to  or  exceeding  the  seismic  magnitude  of  the  explosion 
has  not  yet  been  triggered,  but  because  earthquakes  generate 
wave  types  and  put  energy  in  frequency  bands  that  are  not 
strongly  excited  by  explosions,  there  is  evidence  from  the 
seismic  signals  of  earthquake  generation. 

Geological  Factors  Affecting  Triggering. 

Our  attention  has  been  directed  to  the  question  of  the 
circumstances  under  which  triggering  will  occur  and  the 
environmental  factors  that  determine  the  extent  of  the  fault¬ 
ing. 

Of  the  numerous  questions  raised  by  the  avallaole  field 
observations  of  explosion-induced  faulting,  one  of  the  most 
interesting  is  whether  the  extent  of  faulting  scales  with 
the  yield  of  the  explosion.  The  alternative  is  that  the 
action  is  one  of  triggering  such  that  once  the  yield  exceeds 
some  threshold,  rupture  begins  and  its  extent  is  determined 
by  the  ambient  stress  and  strength  of  the  material  and  is 
independent  of  the  yield.  McKeown  and  Dickey  (1969)  have 
published  data  for  Pahute  Mesa  that  shows  a  straightforward 
relation  of  fault  length  to  yield.  This  dependence,  when 
combined  with  the  absence  of  aftershocks  on  these  faults, 
indicates  that  the  prestress  is  well  below  that  at  which 
slip  would  occur  naturally  and  the  slippage  is  driven  by 


the  3tress  puls**  from  the  explosion.  However,  the  occurrence 
oT  appreciable  slip  and  of  some  aftershocks,  especially  after 
Benham,  is  evidence  that  a  sizable  tectonic  stress  field  is 
present  in  the  rocks  of  Pahute  Mesa . 

The  fact  that  Boxcar,  Benham,  Handley  and  Jorum  all  pro¬ 
duced  visible  surface  faulting  of  about  the  same  extent,  but 
only  Benham  produced  large  numbers  of  aftershocks  has  not 
been  explained.  A  possible  explanation  is  inherent  in  a 
model  of  the  local  structure  that  was  suggested  by  Cummings 
(1966),  Figure  6.1.  He  modelled  the  Timber  Mountain  caldera 
as  a  hole  in  a  homogeneous  elastic  plate  under  uniform  ten¬ 
sion.  He  used  this  model  to  account  for  the  orientation  of 
the  numerous  faults  that  intersect  the  caldera  boundary, 
including  those  in  the  neighborhood  of  the  events  under  dis¬ 
cussion.  Though  the  model  is  undoubtedly  too  simple,  it  does 
account  for  the  orientation  of  the  faults  fairly  well  and  is 
a  credible  representation  of  the  situation. 

In  this  model,  the  large  events  on  Pahute  Mesa  are 
located  roughly  along  the  line  at  right  angles  to  the  direc¬ 
tion  of  regional  tension.  Benham  is  at  about  1.3  caldera 
radii  from  the  center,  and  the  other  shots,  Boxcar,  Jorum, 
Handley,  and  Greeley  are  at  about  2  radii  from  the  center. 
Because  the  effects  of  the  hole  in  a  prestressed  plate  are 
quite  localized,  this  model  calls  for  a  maximum  tensile 
stress  and  a  maximum  shear  stress  at  Benham  about  1.5  times 
those  at  the  more  northern  sites.  Figure  6.2.  The  gradients 
in  stress  are  greater  at  Benham  than  at  the  other  places 


also.  For  each  kilobar  of  regional  tensile  load,  the  maxi¬ 
mum  tensile  stress  is  l.b2  kb  at  Benham  {almost  twice  the 
applied  load)  and  1.22  kb  at  the  northern  site.  The  cor¬ 
responding  values  of  maximum  shear  stress  are  0.73  kb  and 
0.47  kb,  respectively.  The  stress  gradients  at  Benham  are, 
for  the  tensile  stress  0.13  bar/meter  per  kilobar  of  re¬ 
gional  stress,  and  for  the  maximum  shear  stress,  0.07  bar/ 
meter.  At  the  northern  sites  the  corresponding  gradients 
are  about  one-fifth  the  value  at  Benham.  Of  course,  local 
inhomogeneities  may  cause  localized  stress  gradients  over 
short  distances  that  are  much  greater  than  these. 

A  further  point  in  support  of  the  applicability  of  this 
model  is  the  fact  that  the  numerous  Benham  aftershocks  are 
concentrated  along  a  trajectory  whose  direction  is  predicted 
by  Cumming,  for  which  there  was  no  surface  fault,  but  which 
connected  known  segments  of  faults.  Also,  the  late  after¬ 
shocks,  beginning  22  days  after  the  event  are  largely  con¬ 
centrated  to  the  south.  Inside  the  caldera  boundary. 

In  the  absence  of  more  information  about  the  in  situ 
stress  and  the  strength  of  the  faults  it  is  not  possible  to 
Judge  if  a  stress  increase  of  fifty  percent  is  enough  to 
account  for  the  difference  in  response  at  Benham  and  the 
other  shots. 

Further  evidence  of  a  difference  in  Benham  and  Boxcar 
related  to  tectonic  strain  release  is  found  in  the  relative 
excitation  of  SV  and  SH  waves  by  the  two  explosions.  The 
horizontal  component  of  SV  as  observed  at  11  stations  in 
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li¬ 
the  Western  Hemisphere  is  close  to  the  same  for  the  two 
events.  amplitudes,  on  the  other  hand,  are  uniformly 

three  to  four  times  larger  far  the  .enham  event  than  for 
Boxcar,  with  about  the  same  period. 

Milrow  raised  some  interesting  questions,  in  that  it 
produced  so  few  afterevents  in  a  setting  that  is  known  to 
be  highly  seismic.  The  absence  of  aftershocks  later  than 
cavity  collapse  (Engdahl  and  Tarr,  1970)  and  evidence  of 
the  long-term  tectonic  stability  of  Amchitka  Island  (Morris, 
1970;  Anderson,  1970)  Indicate  that  the  rocks  of  that  island 
are  under  less  stress  than  those  of  the  Basin  and  Range 
Province.  One  Inference  is  that  the  island  is  mechanically 
decoupled  from  the  active  seismic  region  in  which  it  is 
embedded.  Anderson  interprets  "the  deformation  of  the  sur- 
flcial  rocks  of  the  central  Aleutian  Ridge  as  related  to 
volcanlsm  and  plutonism,  and  as  such  the  deformation  does 
not  reflect  in  a  first  order  sense  the  regional  compresslonal 
stress  across  the  arc." 

The  strain  energy  released  in  the  Milrow  aftershocks 
almost  certainly  came  from  the  static  strain  field  produced 
by  the  explosion  and  subsequently  relieved  by  the  collapse. 

Thus,  the  prompt  release  of  tectonic  strain  by  an  ex¬ 
plosion  is  governed  by  a  number  of  factors  that  can  be 
evaluated  in  the  field  but  have  not  been  adequately  treated 
yet.  These  include  the  magnitude  and  orientation  of  the 
regional  stresses  and  the  stress  concentrations  in  the 
neighborhood  of  the  shot  caused  by  the  local  geology.  The 
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occurrence  of  an  unstable  rupture  Hint  would  give  rise  to  a 
large  earthquake  seems  to  require  that  the  explosion  take 
place  near  a  fault  that  is  stressed  almost  to  the  point  of 
failure.  In  this  event  the  identification  ol'  the  signal  as 
coming  from  an  explosion  might  ee  very  difficult. 

It  has  been  recognized  that  surfs  e  wave  magnitudes 
from  NTS  events  are  systematically  Isrgerl  relative  to  mb, 
than  those  for  explosions  In  Central  /e-la,  Novaya  Zemlya 
and  the  Sahara.  For  example,  for  a  large  Novaya  Zemlya 
event,  mb  =  6.3,  the  ratio  m^  :  is  I.30.  Tne  average 

for  two  Central  Asia  events,  mL  -  o.O  and  3. hi  the  ratio 
is  1.32.  For  one  Sahara  event,  m^  =  p.7,  the  ratio  is  1.32. 
For  events  Greeley  and  Boxcar  the  ratios  are  both  1.16. 

0.  Nuttli  has  determined  shear  wave  magnitudes  and  found 
for  these  events  that  the  shear  body  wave  magnitude  is 
close  to  Mg  for  all  sites,  the  ratios  '«Shear  :  Ms  ran^ihg 
from  1.02  to  I.I3.  The  relatively  high  values  for  both 
M3  and  for  NTS  events  have  been  Interpreted  as  due 

to  the  contribution  of  tectonic  strain  release.  The  con¬ 
tribution  to  the  shear  wave  magnitude  from  the  SH  component 
is  very  small  for  most  of  the  events  examined.  A  detailed 
study  of  the  Benham  shear  wave  magnitude  is  still  to  be 
done . 

Aftershock  Excitation. 

In  collaboration  with  J.  T.  Cherry,  Jr.,  of  Livermore 
Radiation  Laboratory,  Livermore,  California,  a  mechanism  by 
which  an  underground  explosion  may  trigger  a  swarm  of  small 


earthquakes  has  been  Investigated.  The  effect  to  be  explained 
Is  the  permanent  alteration  of  the  tectonic  environment  in 
the  zone  usually  considered  to  respond  elastically  to  the  ex¬ 
plosive  loading.  A  non-elastic  effect  must  be  produced  be¬ 
yond  the  usually  defined  zone  of  cracking,  crushing,  etc., 
one  that  causes  small  earthquakes  to  occur  over  an  extended 
period  of  time. 

The  suggested  mechanism  Is  the  interaction  of  elastic 
shear  waves  proceeding  outward  from  the  region  around  the 
explosion  with  localized  gradients  in  the  ambient  tectonic 
stress  field.  As  a  result  of  the  interaction  a  field  of 
small  slipped  areas  or  dislocation  loops  is  formed  and  the 
continued  action  of  the  ambient  field  on  these  dislocations 
results  in  the  swarm  of  earthouaket .  The  shear  waves  may 
proceed  directly  from  the  explosion  or  be  produced  by  re¬ 
flection  of  compress ional  waves  at  the  free  surface  or  other 
boundaries.  The  stress  concentrations  are  the  result  of 
either  material  inhomogeneities  or  the  effect  of  preexisting 
cracks  in  the  medium,  which  is  in  equilibrium  in  an  ambient 
regional  stress  field. 

Our  appeal  to  the  role  of  local  stress  concentrations 
follows  closely  the  hypotheses  of  Droste  and  Teisseyre 
(1959,  i960)  and  of  Mogl  ( 1962 )  regarding  the  preparation 
of  a  region  for  earthquakes  and  aftershocks. 

The  equations  of  motion  for  an  elastic  medium  under 
initial  stress  are  derived  by  Biot  in  The  Mechanics  of  Incre¬ 


mental  Deformation.  If  the  body  forces  and  initial  stress 
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are  all  zero,  we  have  the  usual  equations  of  motion.  For 
the  body  forces  zero  and  the  Initial  stresses  constant.  Blot 
has  shown  that  P  and  S  waves  can  propagate  In  the  medium 
and  that  the  velocities  depend  on  the  amount  of  the  Initial 
stress, and  an  Initially  Isotropic  medium  is  in  general  not 
Isotropic  under  the  ambient  stresses. 

If  the  body  forces  are  zero  and  if  the  coordinate  axes 
are  taken  as  the  principal  axes  of  the  ambient  stress,  but 
the  ambient  stress  Intensity  is  permitted  to  vary,  the  two 
equations  are  coupled  through  the  shear  strain  term,  which 
Involves  the  stress  gradients  as  coefficients.  Although  the 
solution  of  the  resulting  equations  Is  very  difficult.  It  can 
be  easily  seen  by  substitution  that  a  pure  plane  shear  wave 
propagating  parallel  to  one  of  the  coordinate  axes  is  not  a 
possible  solution.  For  example,  a  plane  shear  wave  travelling 
parallel  to  the  y  axis,  displacement  parallel  to  the  x  axis, 
must  be  accompanied  by  a  particle  acceleration  parallel  to 
the  y  axis.  A  pure  congressional  wave,  on  th°  other  hand, 
can  exist! 

Physically,  the  stress  gradient  acts  analogously  to  the 
tension  in  a  stretched  string,  resulting  in  an  equivalent 
spring  that  converts  a  displacement  into  an  acceleration. 
Because  the  stress  gradient  is  a  multiplier  of  the  strain  in 
the  equation,  a  small  strain  can  produce  a  large  acceleration 
if  the  stress  gradient  is  large,  and  a  change  in  sign  of  the 
stress  gradient  will  change  the  sign  of  the  acceleration. 

Large  stress  gradients  are  easily  produced  over  the  short 
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We  propo3«*  tnen  that  a3  a  shear  wave  passes  through  a 
zone  In  whl  -n  a  strong  stress  gradient  changes  sign,  a  re¬ 
sultant  tensile  stress  may  act  momentarily , reducing  the 
eiTe 'tlve  pressure,  and  slip  may  occur  under  the  action  of 
the  ambient  stress.  The  slipped  areas  created  in  this  may 
become  the  nuclei  for  growth  of  the  dislocations  under  the 
action  of  the  ambient  field,  in  the  manner  proposed  by  Droste 
and  Teisseyre. 

The  Benham  aftershock  swarm  exhibits  features  that  may 
be  explained  by  this  mechanism.  Pig.  6.3  shows  the  epi¬ 
centers  and  the  hypocenters,  on  a  vertical  radial  plane,  of 
the  first  19  aftershocks  reported  by  Hamilton  and  Healy 
(I9b9).  These  are  only  a  small  sample  of  the  many  events 
during  the  first  six  hours  after  Benham,  but  they  show  that 
early  aftershocks  occurred  at  three  to  six  kilometers  from 
the  shot,  and,  in  fact,  outline  the  zone  of  most  of  the  ac¬ 
tivity  for  the  first  three  weeks  after  the  event.  The  mechan 
isn  by  which  the  region  wa3  primed  for  the  aftershocks  acted 
quickly  and  did  not  move  out  slowly  from  the  explosion  site. 

Pig.  6.4  is  a  hypocenter  plot  of  all  668  events  located 
by  the  USGS  through  April  18,  1969,  on  a  single  plane.  Super 
Imposed  are  two  rays:  one  defines  the  deepest  ray  for  a  P  to 
S  conversion  by  surface  reflection  for  which  more  than  half 
the  reflected  energy  is  in  the  form  of  S  waves  for  an  ideal¬ 
ized  medium.  All  rays  above  this  one  out  to  large  angles  of 
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Incidence  at  the  surface  exhibit  t  .is  property.  The  ray 
appears  to  form  a  rough  lower  envelope  to  the  activity.  The 
velocity  model  used  by  the  USGS  was  used  in  drawing  this  ray. 

The  other  ray  represents  an  SV  wave  taking  effect  about 
is  to  the  vertical  at  the  source.  This  is  the  angle  at 
which  a  beam  of  energetic  S  waves  should  radiate  from  a  hori¬ 
zontal  depth  force  without  moment  and  was  selected  on  the 
basis  of  a  model  of  a  near-surface  explosion  proposed  by 
Gupta  and  Kissllnger  in  1966.  This  ray  also  passes  along 
the  lower  bound  of  the  activity.  For  either  source  of  S- 
waves,  the  cut-off  of  activity  below  about  5.5  km  may  then 
be  an  effect  of  the  sharp  velocity  contrast  at  5  km  depth. 

The  circled  points  are  hypocenters  of  late  events  in  the 
swarm,  mostly  south  of  the  regional  swarm  regions. 

The  distribution  of  the  aftershocks  shown  by  Hamilton 
and  Healy  (1969)  raises  the  question  of  why  the  events  are 
distributed  as  they  are.  The  absence  of  activity  to  the 
east  suggests,  under  the  present  hypothesis,  that  the  main 
fault,  which  did  slip  at  shot  time,  effectively  decoupled 
shear  waves  from  the  eastern  zone.  The  concentration  of 
events  between  previously  mapped  fault  segments  suggests 
that  the  swarm  has  acted  to  complete  the  formation  of  a  linear 
subsurface  fracture.  If  so,  the  swarm  has  some  similarities 
to  natural  earthquake  swarms,  for  example  that  in  Matsushiro, 
Japan.  There  a  great  swarm  accompanied  the  formation  of  a 
new  fault,  according  to  Japanese  geologists,  activity  died 
out  in  the  original  active  zor.c  as  it  did  in  Benham,  and 
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In  later  stages  of  the  swarm  the  activity  expanded  outward 
Into  fresh  material. 

The  limited  evidence  indicates  that  moderate  to  strong 
earthquakes  occur  on  preexisting  faults,  and  no  case  is  known 
to  the  author  in  which  a  great  earthquake  was  accompanied  by 
the  formation  of  a  new  fault  in  virgin  material.  It  may  be 
that  the  formation  of  a  new  fault  is  normally  accompanied  by 
an  earthquake  swarm,  perhaps  initially  microearthquakes, 
though  the  converse  is  by  no  means  indicated,  that  is,  that 
all  swarms  are  associated  with  fault  formation. 
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Figure  6. 


Figure  6. 


Figure  6.3 


Figure  6.4 


Figure  Captions 

Model  of  Pahute  Mesa  -  Timber  Mountain  caldera 
tectonics  after  D.  Gumming,  cr"  is  the  re¬ 
gional  tensile  stress.  Nuclear  events  are 
designated:  H  -  Handley,  J  -  Jorum,  G  - 

Greeley,  Bo  -  Boxcar,  Be  -  Benham.  The  dashed 
curves  are  predicted  orientations  of  faults. 
Variation  of  maximum  tensile  stress, ^7  ,  and 
maximum  shear  stress,  S,  with  distance  to  the 
north  from  the  caldera  boundary.  Distance  is 
in  terms  of  caldera  radii. 

Epicenters  (left)  and  hypocenters  of  the  first 
19  Benham  aftershocks  listed  by  the  USGS.  Body 
wave  magnitudes  are  in  parentheses. 

Hypocentral  plot  of  Benham  aftershocks  through 
April  18,  1969.  The  solid  ray  is  a  P  to  S 
conversion  at  the  smallest  angle  of  incidence 
for  which  more  than  one-half  the  reflected 
energy  is  in  the  shear  wave.  The  dashed  ray 
represents  the  direction  of  the  maximum  ampli¬ 
tude  of  SV  waves  generated  by  the  horizontal 
component  of  the  explosive  loading  for  a  near¬ 


surface  source. 


adapted  from  Cummings  (1968) 


10  km 


7.  P  Pulses  from  Earthquakes  Triggered  by  Explosions. 

Rene  Rodriguez 


Introduction 

The  purpose  of  this  research  was  to  determine  in  the  time 
domain  the  influence  of  an  earthquake  triggered  by  explosions 
into  the  final  P  pulse  shape.  Two  kinds  of  sources  were  con¬ 
sidered  to  act  in  a  region,  namely:  an  explosive  source,  and 
a  tectonic  source  (moving  fault)  triggered  by  the  arrival  of 
a  signal  from  the  explosion. 

a)  Explosive  Source.  An  explosive  point  source  in  an 
homogeneous,  isotropic  media  can  be  represented  by  three  mu¬ 
tual  perpendicular  couples  of  equal  point  forces  without 
moment  (Fuchs,  1966).  The  total  effect  can  be  represented 
by  a  superposition  of  the  three  couples,  because  of  the  linear 
dependence  between  displacements  and  forces.  This  is  equiva¬ 
lent  to  considering  the  explosion  and  its  effects  only  in  the 
elastic  zone,  where  the  infinitesimal  strain  theory  is  valid 
(Kisslinger,  1963).  From  Nakano  (1923)*  the  displacements 
due  to  an  explosive  point  source  were  calculated  and  they  are: 
in  cartesian  coordinates: 


Up  •  JUL—  1  (  **♦  jj**«.i»x)i  4.  (t*y  j  fi.)  C1) 

41T  f  l  J 

in  polar  coordinates: 


±ir  a 


K  fc 

4-TTf  woy? 


fir 


(2) 
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where  compress Iona 1  wave  velocity 

R  =  distance  between  source  and  receiver 
f  «  density  of  the  media 

<€  -  equivalent  moment  of  the  couple  (strength  of  the 
explosive  source) 

§r  =  radial  unit  vector. 

b)  Tectonic  Source.  A  .oving  double  couple  is  utilized 
to  model  an  extending  source.  Following  Berckhemer  (1966), 
consider  the  fault  plane  in  the  xz  plane  and  a  double  couple 
acting  at  the  origin  0  (Fig .  7.1).  Let  B  be  a  point  on  the 
focal  plane  given  by  it  cartesian  coordinates.  Let  the  force 
per  unit  area  K  be  acting  over  the  surface  element  df  , 
then  the  resultant  force  will  be  given  by 

<iK  *  JU!  elf  (4) 

Let  the  receiver  A  be  at  a  distance  r  from  the  source,  very 
large  compared  with  the  diameter  2r'  of  the  focal  area.  The 
rays  from  all  the  source  elements  to  A  will  nearly  be  parallel. 
This  approximation  can  be  Justified  if 

f*  <<  Amin  r 

where  is  the  wave  length. 

If  the  focal  process  starts  at  the  origin  0,  the  stress  re¬ 
lease  at  any  other  point  such  as  B  will  start  after  a  delay 
time 

(5) 


where 


»»  time  due  to  the  fracture  process. 


The  delay  time  at  the  receiver  A  will  be: 


U  .  t(Cs)  -  r 

where: 

fir  -  JL 
r 

$  »  fi  ♦  r  k 

ti  «  _5_  .  Co*  G 

It  i  ,  5m© 

I* »  «  -i-  s  5m  8  t0*  ^ 
r  r 

r  .  i*)'4 


(6) 


®  ^  for  P  waves  v>%  for  S  waves 

In  each  source  element  a  double  couple  acts.  Integration 
over  the  total  focal  plane  leads  to  the  resulting  pulse  of 
the  displacement  at  A.  Prom  the  expressions  given  by  Nakano 
(1923)  the  displacement  at  A  will  be 

it'  ■  “lyfej+jaL*  H(rt)  dta;  (9) 

u«  .  f*f‘  a,  aj  (10) 

4tff  r  of  Je  'd 

/ 

where  Up  stands  for  P  waves,  Us  for  S  waves,  a,  b,  c,  d,  are 
the  dimensions  of  the  fault  and  is  the  tangential  unit 
vector. 

The  composite  source. 

Consider  an  isotropic,  homogeneous  media.  Let  the  explo 
slon  be  detonated  at  a  certain  time  t^.  Waves  spreading 
radially  will  travel  around  the  explosive  source  in  all  direc 
tlons.  It  is  assumed  as  a  working  model  that  as  soon  as  the 
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first  wave  front  from  the  explosion  reaches  the  fault  area 
in  one  point,  this  will  be  the  starting  point  for  the  trigger¬ 
ing;  this  will  occur  after  a  delay  time  (Pig.  7.2) 

t  -  -2.  (id 

«r 

where  M  is  the  distance  between  the  explosion  and  fault. 

Any  other  point  in  the  fault  such  as  B  will  start  releasing 
energy  after  a  time  delay 

l-  (12) 

At  this  point  if  we  consider  the  fracture  process  to  be  linear 

(unilateral  or  bilateral)  propagating  in  +  x  direction 

-|-  (13) 

where  C  -  fracture  velocity. 

Then  the  total  time  delay  from  the  explosion  to  any  point  on 
the  fault  will  be: 

<»> 

Meanwhile  direct  waves  from  the  explosion  traveling  thru  the 
medium  will  reach  the  receiver  after  a  time  delay: 


and  waves  from  the  fault  after  a  time: 


t-  (fr.f  (15) 
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when  the  fracture  process  is  unilateral  or 

1  (  *  .  T  .  r  +  (?€*• » iSmBtouQ  \  fj 

t  "  V  1?  T  3  f 


(16) 


when  it  is  bilateral. 

V  •  for  P  waves  V  =  Vs  for  S  waves 

Then  the  total  displacement  due  to  the  composite  source  at 

the  receiver  t-  will  be 


Vf 


-«■*»  . 

♦IT  f  R  yJf 


(17) 


u* 


*C 


( w*.w.y Thant  (is) 


If  we  assume  further  that  the  radius  of  the  composite  source 
r"  is  very  small  compared  with  the  distance,  which  can  be 
satisfied  if 


r*  <<  Xw«*r 


(19) 


then  the  displacements  can  be  written  as: 


r  JLfc-, —  F (i-  B.)  4  .si  C  (b  clT  dr 

41TJR^  11  V  tr  )6  )4  ' 


(20) 


u»8  Jfi 


r</r‘1  dfdt  (21) 

4ir  f  r  o{  Je 


for  unilateral  fracture,  and 


Up,  J<lL»f(1.4)  +  JUAaSatZukC  i(Xl)  dtdt 

♦irfRo*  V  iirfrdj 


(22) 
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3 


Kfi 


i *■"»>  * ( i.  w)}Vt  f*  f  *  ; (X\ }  diar  (23) 

♦tfroj  J*  Jd 


for  bilateral  fracture. 

For  the  final  evaluation  of  the  displacements  an  assumption 
about  the  time  function  for  P  and  S  wave  forms  at  the  source 
has  to  be  made . 

1.  Time  function  for  the  explosive  source. 

Two  kinds  of  functions  will  be  assumed: 

Type  1  Explosion  (EX).  The  displacements  at  large  distance 
corresponding  to  step  in  pressure  at  the  source  is  given  by 
(Kisslinger,  1963) 

u  ,  u,  e**T  ta.r 


which  is  a  damped  sinusoidal  motion  propagating  radially  with 
angular  frequency  wQ  where: 


Assuming  for 
and  for  Vp  - 


t . 

t-  u 

(24) 

o<  = 

m  Wo 
( m  •»  m*)’* 

(25) 

U)|; 

*0f 

(26) 

a,  (i+a m) 

ai=  radius  of  the  cavity 

JL  «  ratio  of  shear  modulus  to  Lame's  constant. 

x 

the  earth  that  Poisson's  ratio  is  0.25,  m  =  1 
7  km/sec 


w»-  i±JL _ L_ 

3  <>i 


Q<  « 


14  J_ 
3  <*i 
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For1  a  ve^y  large  cavity  radius,  for  example  a  »  1  km 

u3°  =  6.6  1/sec 
=4.66  l/sec 
fr«<?  =  1.05  Hz 

The  corresponding  normalized  displacement  for  this  type  of 
time  function  is  plotted  in  Figure  7.3  taking  of  as  a  parameter. 
Type  2  Explosion  (EX).  A  step-like  function  at  the  source 
will  be  assumed  as  a  second  type  of  time  function  for  an  ex¬ 
plosion  of  the  form: 

u  ,  ut  t  e-*T  (27) 


where  will  be  evaluated  by  comparison  during  the  numerical 
evaluation  of  the  pulses.  This  type  of  function  is  similar 
to  the  pressure  function  calculated  by  Toksoz  (1967)  from  sur 
face  wave  studies  of  the  Benham  nuclear  explosion.  The 
value  of  1.5  for  in  the  time  function  seems  to  fit  best  the 
observed  Love  wave  spectrum. 

Plots  of  the  normalized  displacement  for  different  alpha's 
are  shown  in  Figure  7.4. 

2.  Time  function  for  the  tectonic  source. 

Various  kinds  of  wave  forms  have  been  proposed  by  many  dif¬ 
ferent  authors.  Bollinger  (1967)  gives  a  list  of  the  P  wave 
forms  used  most  often: 


Boxcar 


Step  function 


Damped  sinusoid 


0 


Uft) 
e**1*  Sin  (If 


t  >0 
t  <0 


(28) 


H  U  (-o 


(*t 


Prom  spectral  studies  of  body  waves  by  Ben-Menahem  et  al 
(1965)1  the  gross  structure  of  the  source  time  function  was 
found  to  be  of  the  form: 

( 1  -  e-*x) 

with  o  c  1  i  10  (29) 

(* 

Studies  on  the  observed  surface  waves  by  Ben-Menahem  and 
Toksoz  {1963  a,b)  and  Stauder  and  Hirasawa  (1965),  show 
that  the  time  function  at  the  source  can  be  expressed  by  a 
step-like  function.  According  to  their  results  the  follow¬ 
ing  particular  form  can  be  assumed  for  the  time  function: 


1-  e-* 

t*0 

(30) 

,  0 

t<  0 

or  in  dimensionless  time  form 

f 

1 .  e-QT 

T  *  O 

K  T)  . 

(31) 

1 

0 

T  <  O 

where 


T  s  ft 

and  b  (32) 

Q  =  -*£- 

b  =  fault  length 
Vs  «*  S  wave  velocity 


Wave  forms  of  P  waves  are  plotted  considering  Q  as  a  param¬ 
eter  (Pig.  7.5).  Ben-Menahem  and  Toksoz  presented  the  value 
of  (3*  1/22  from  surface  waves  of  the  1952  Kamchatka  earth¬ 
quake.  Stauder  and  Hirasawa  present  a  corresponding  value 
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for  Q  =  10.  Prom  the  considerations  of  Stauder  and  Hirasawa, 
Q  may  be  considered  to  be  almost  constant  for  any  earthquake 
rather  than  constant,  because  it  seems  to  take  a  longer 
time  to  release  the  stored  strain  energy  in  a  larger  source, 
than  in  a  smaller  one. 

In  the  present  paper  two  kinds  of  time  functions  will  be 
assumed  for  P  waves: 

Type  1  Fault  (FAU).  A  damped  sinusoid. 

{  (t)  =  e-<*T  5m(3T  (33) 


will  be  assumed  for  P  wave  forms  from  earthquakes,  where  the 
value  of  (3  will  be  chosen  by  comparison  during  the  numerical 
computations  of  seismic  pulses. 

Type  2  Fault  (FAU).  A  step-like  function  for  the  source 
function  as  given  by  Ben-Menahem  (1963  a,b)  and  Hirasawa  and 
Stauder  (1965) 

f  (x\ 

l  h  #  4  M 

(34) 


t-e* 

0 


T  »o 
t  <0 


The  reason  for  assuming  this  time  function  instead  of  a  more 
complete  one,  for  example  that  given  by  Ben-Menahem,  et  al 
(1965),  is  because  we  want  to  gain  simplicity  for  the  evalua¬ 
tion  of  the  integral  over  the  fault  area,  which  for  the  case 
of  a  function  involving  Heaviside  function  cannot  be  evaluated 
analytically  as  we  wish. 
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Models . 

Theoretical  studies  by  Berckhemer  and  Jacob  (1965) ,  Savage 
(1965*  1966),  Stauder  and  Hlrasawa  (1965)  show  that  a  uni¬ 
lateral  fracture  process  has  more  influence  on  seismic 
pulses  than  bilateral  faulting.  By  the  nature  of  the  prob¬ 
lem,  unilateral  faulting  seems  the  more  suitable  process  to 
model  a  triggered  earthquake  although  bilateral  faulting  is 
also  possible.  Because  we  are  interested  in  seeing  the  in¬ 
fluence  of  triggered  earthquakes  on  seismic  pulses  from 
nuclear  explosions,  a  unilateral  fracture  process  will  be 
assumed  in  the  development  of  the  final  displacement  formulas. 
Model  1.  Type  2EX  +  Type  3.PAU  (Fig.  7.6) 

f(t)  x  e- uZ 

!  (X)  z  e**T  S.nf»t  (35) 

With  these  two  time  functions,  expressions  for  the  displace¬ 
ments  corresponding  to  formulas  (20)  and  (21),  after  inte¬ 
gration  over  the  fault  area,  becomes: 


UP  = 


o4  ti  e-  Ml1 

oa  | .  e- l  ♦  co«  (rti]  *  e-  ♦  co* |ti»]  + 

(36) 


u*  =  o»  j.  e-^r  <•*  * 

♦  1  *•«»  -  e  ^1*  J 


(37) 
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Where 

0  -  Kl£t  q9.  K(d»S  n  ■  ^ 

*  4trrr^  <*  (38) 

tj.  1.  a.  t*. 

1  <jf  *®P  <H  ' 


T>‘  »-(-$♦■£ -^-i^) 


t?  = 


t- 


M  _  t 

®p"  ®1 


tJs 

X4  s 
tj. 


i-  (*-$ 

1  -  ( %  *  Vb  ^ ^  -  2S2$fttJt) 


t.=  14  -b(l--S^) 

T9  =  tT  -  b(l  -  £B|J) 

(39) 


V  >  Vp  for  P  waves  and  V  ■  V8  for  S  waves. 


Model  2.  Type  1EX  +  1FAU  (Pig.  7.6) 

F(TW  e‘°rt  Smufct 

f(Tl  e*^T  Sir> fit 


m 


With  the  same  configuration  as  Model  1  the  integrated  ex¬ 
pressions  corresponding  to  formulas  (20)  and  (21)  with  the 
time  functions  assumed  are: 


Up-  Oj  6*  ^  0|  |  “  ^  ^  [fc*(ltt+C«|(ltil  4. 

+  e-*T>  (*.«*»*$ 

.  g-  [  4iH(»T»  4  Cm  jrt*]  | 


(41) 


u*  ,  0»  I*  C"  Cm^T*]  ♦  0’^’  [*t«<itT4  Ukfftrl  4 

4  C  %mp>‘T|  4 Cm |*t§]  .  6^^  ^ 


(42) 
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Where  Dl,  D2  and  D3  are  the  same  as  In  Model  1. 

Model  3.  Type  1EX  +  Type  2FAU  (Pig*  7*6) 

FLt)  z  e-  0,t  (43) 

i  &  ri- 

Using  again  the  same  configuration  of  Model  1  the  integrated 
expressions  corresponding  to  formulas  (20)  and  (21)  with  the 
time  functions  assumed  are: 

Uf  =  0,  e-*“  *  D,  j  .  e-**»  -  e-pU  *  e'"T‘  }  (44) 

Ul!  0,  1  .  e-*'1’  -  *  e-11''*]  (45) 

are  given  by  formula (39) 

*1 .  r*  1 .  • 

Where  again  Dl,  D2,  and  D3  are  the  same  as  in  Model  1. 

It  is  noticed  that  whatever  the  time  function  is,  we  have 
the  same  time  delays,  namely 


T‘*  -*> 

T1  t  -fl.  *  S.  .  9.Smj8  j 

*  JEL ♦  X- 

*  ^ 

T4  *  -  *-*1  *(m  4 


(46) 


for  P  waves,  and 


140 


-pfc  .  m  t  r  _  a  Sin 6  Cot 

Dp 


Tj  ^  f±  +  JL_ 

<5p  dj 


T  8  =  I2_  +  JL  4-  b<  1  -  SftiS  >  -  Q  ^in9  Coftf 

*p  C*  e.  at  '  vJfc 


T»  ■  + 


for  S  waves. 

The  sense  of  the  first  (direct)  motion  from  the  explosion 
will  always  be  up  if  the  positive  sign  is  taken  in  the  direc¬ 
tion  of  Z  axis.  The  second  arrival  will  correspond  to  the  P 
wave  triggered  by  the  explosion  due  to  the  contribution  of  T2 
or  T3,  whichever  is  minimum.  01  course  if 

a  Srnd  C6&0  <1  T2  is  minimum,  or  if 


a  Sine  coipf  >  i  T3  is  minimum,  depending  in 

either  case  on  the  values  of  0  and  ^  .  This  minimum  value 
also  will  determine  the  polarity  of  the  second  arrival  for  a 
given  constant  distance,  because  the  difference  in  sign  be¬ 
tween  T2  and  T3.  The  difference  between  T2  and  T3  is  strong¬ 
ly  dependent  on  the  width  of  the  fault  considering  that  •  ,  , 

and  Vp  are  fixed,  that  is: 


n*i  -  iTji 


(47) 


On  the  other  hand,  because  the  difference  in  sign  a  reverse 
in  polarity  will  occur  after  the  contribution  of  the  third 
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term  (one-fourth  of  a  complete  cycle),  therefore  the  theory 
predicts  a  possible  way  of  estimating  the  width  of  the  fault 
independent  of  the  assumption  for  the  time  function,  if, 
of  course,  in  practice  this  second  arrival  could  be  identi¬ 
fied  with  accuracy  from  actual  records. 

A  third  spike  is  expected  due  to  the  contribution  of  T4  or 
T5,  whichever  is  minimum,  where  the  length  of  the  fault  plays 
an  important  role.  This  arrival  will  correspond  to  the 
stopping  phase  (Savage  1965,  1966).  The  same  analysis  can 
be  made  for  the  S  wave  generated  at  the  fault. 

Numerical  Calculations. 

Fortran  programs  have  been  written  to  evaluate  the  displace¬ 
ments  at  the  receiver  for  all  three  models  and  variable  param¬ 
eters  such  as  distance,  relative  orientation  of  the  sources, 
relative  strength,  etc. 

Model  1. 

These  programs  calculate  the  P  and  S  pulses  at  the  receiver 
due  to  the  composite  source  for  different  parameters,  namely: 

a)  r  =  distance  between  fault-explosion  to  station 

b)  Q  and  i  related  the  orientation  of  the  station  with 

respect  to  the  fault 

c)  M  *  distance  between  explosion  and  fault 

d)  and  related  the  orientation  of  the  explosion 

with  respect  to  the  fault 

e)  b  »  length  of  the  fault. 

At  the  beginning  values  of  °<  *  0.5  and  (3  -  0.1  were  assumed 
arbitrarily,  considering  only  that  the  rise  and  decay  for  the 
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explosive  source  should  be  faster  than  for  the  propagating 
fault.  All  the  early  results  for  Model  1  correspond  to 
these  two  values.  It  then  became  clear  that  these  coeffi¬ 
cients  were  critical,  and  various  attempts  to  calculate  the 
order  of  these  have  been  made.  Attempts  to  relate  these 
coefficients  to  the  magnitude  and  distance  and  the  time  of 
duration  of  observed  pulses  were  unsuccessful.  More  suc¬ 
cessful  was  an  attempt  to  relate  with  (*  directly,  equat¬ 
ing  the  theoretical  displacements  due  to  an  explosive  point 
source  acting  alone  and  the  displacement  due  to  a  double 
couple  of  equivalent  strength.  The  numerical  values  obtained 
have  been  used  to  calculate  the  values  as  shown  in  Pig.  7.7. 
As  predicted  by  the  theory,  the  first  arrival  corresponds  to 
the  direct  pulse  from  the  explosion,  the  second  to  the 
triggered  earthquake,  the  third  is  the  stopping  phase  from 
the  earthquake,  and  the  fourth  the  S  arrival.  The  differ¬ 
ence  in  frequency  content  between  the  pulse  from  the  explo¬ 
sion  alone  and  that  corresponding  to  the  fault  is  striking. 

In  these  calculations  the  following  values  for  the  different 
parameters  were  taken: 
o(  «  5.0 
f*  -  0.7 

o'? 

Strength  of  the  tectonic  source  “.s.w,.  io  dyn-cm 
Distance  fault-station  r  *  100  km 

Distance  explosion-fault  M  «  15  km 

Fault  length  b  *  10  km 

Strength  of  the  explosive  source  10^  dyn-cm 
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e  -  30° 

4  =  120 

©i-  150° 

4i-  30' 

Pig.  7.8  is  the  same  as  Pig.  7.7  but  with  strength  of  the 
explosive  source  decreased  to  10  dyn-cm. 

Model  2. 

The  program  for  this  model  calculates  the  displacement  due 
to  a  composite  source  as  assumed  in  Model  2.  Displacements 
have  been  calculated  ut;ing  the  expressions  derived  by  Kiss- 
linger  (1963)  and  by  amplitude  equalization.  As  in  Model  1 
we  have  the  four  spikes  and  more  or  less  equal  frequency  con¬ 
tent  .  In  this  case  also  the  same  parameters  of  Model  1  have 
been  used  (Pig.  7.9). 

Model  3. 

This  program  calculates  the  displacement  due  to  a  composite 
source  as  assumed  in  Model  3.  The  value  of  ,  as  suggested 
by  the  theoretical  and  observational  calculations  of  Kiss- 
linger,  was  taken  as  4.66,  which  corresponds  to  a  very  large 
cavity  radius  (1  km) .  Beta  again  was  calculated  by  amplitude 
equalization.  In  this  model  the  arrivals  observed  in  Models  1 
and  2  are  clearly  seen.  The  frequency  content  is  slightly 
different  from  those  models,  probably  because  of  different 
values  of  the  coefficients  used.  This  model,  not  only 
theoretically  but  also  experimentally  seems  to  fit  better 
the  actual  physical  problem  (Pig.  7.10). 

Currently  we  are  trying  to  improve  the  Portran  programming 
to  obtain  separately  first  pulses  only  from  the  explosion  and 
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from  the  triggered  earthquake,  as  well  as  the  combination 
of  both.  This  will  give  a  better  picture  of  what  we  are 
intending  to  do.  After  this,  the  next  step  will  be  to 
analyze  the  relative  importance  of  each  one  of  the  param¬ 
eters  used  in  the  theory  into  the  final  p  pulses.  Further 
on  all  the  theoretical  P  pulses  calculated  will  be  con¬ 
volved  with  the  crust-seismograph  response  to  determine  the 
actual  theoretical  seismogram.  In  this  sense  a  convolution 
Fortran  program  has  been  already  written  and  the  response 
for  the  WWSSN  seismograph  and  crust  response  has  been  already 
obtained  in  digital  form,  taking  the  models  given  by  Ben- 
Menahem. 

Also,  some  practical  applications  of  this  theory  will  be 
applied.  For  this  purpose,  P  waves  from  the  Boxcar  nuclear 
explosion  are  being  analyzed,  where  very  peculiar  kinks  in 
the  first  cycle  have  been  found  at  some  of  the  standard  sta¬ 
tions.  What  those  kinks  mean  physically  we  are  not  in  a 
position  to  say  yet;  however,  on  the  basis  of  a  better  under¬ 
standing  of  the  theory,  perhaps  we  will  be  able  to  answer 
this  very  interesting  question. 

The  analysis  developed  here  can  easily  be  extended  to 
the  case  in  which  the  triggering  occurs  at  the  time  of 
arrival  of  an  S-wave  proceeding  from  the  explosion.  The 
details  of  this  extension  will  depend  on  the  model  chosen 
for  the  generation  of  these  S-waves,  e.g.  by  reflection,  by 
crack  formation,  etc. 
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8. 

A  study  of  smaller  BENHAM  aftershocks,  begun  on  a 
project  with  the  U.S.  Geological  Survey  during  the  summer 
of  1969>  was  completed  under  the  Contract  during  this  report 
period,  and  a  paper  prepared  for  publication.  The  title  of 
the  paper  is  "Smaller  Aftershocks  of  the  BENHAM  Nuclear  Ex¬ 
plosion,  "  by  William  Stauder.  The  Abstract  follows: 

The  examination  of  the  BENHAM  aftershocks  is  extended 
to  162  earthquakes  approximately  an  order  of  magnitude  smaller 
than  those  previously  studied.  The  spatial  distribution  of 
these  smaller  events  is  similar  to  that  found  for  larger 
events  examined  by  Hamilton  and  Healy.  Focal  mechanisms  are 
also  similar  to  those  of  the  larger  aftershocks:  dip-slip 
along  northeasterly  trending  zones,  generally  strike-slip 
along  a  north-south  trending  zone  west  of  the  shot-point. 

About  one-third  of  the  mechanism  solutions  are  ambiguous, 
capable  of  either  interpretation.  Amplitude  data  are  suc¬ 
cessfully  used  to  resolve  the  ambiguity  in  about  half  these 
cases  by  selecting  the  solution  which  gives  a  notably  lesser 
value  of  the  variance,  ^(A^  Qt)S  -  kA.^  caic)2/N  of  the  amP11" 
tude  residuals.  The  value  of  k  is  related  to  the  magnitude 
of  the  shocks.  A  b  value  of  1.09  indicates  that  aftershocks 
of  the  explosion  follow  a  recurrence  rate  normal  for  regional 
earthquakes  or  aftershock  sequences. 


Focal  Machanloaa 
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9.  An  Application  of  P  Hm  Stationary  Phase  Approximation  In 
Determining  Source  Parana  tart  by  equalisation  Procedures. 

W.  Staudar  and  M  U 11  Una 

A  atlaulua  for  tha  uaa  of  P  Mva  aaplltuda  data  In 
determining  aourca  parameters  arlaaa  fro*  tha  paraonal  obser¬ 
vation  of  similarity  of  wave-form  of  tha  P  wave  fro*  station 
to  station  for  a  given  aarthquaka  and  of  tha  diminishing  of 
aaplltuda  of  tha  inltUl  P  wave  at  stations  which  after 
analysis  and  ooaparUon  of  all  tha  data  ara  shown  to  ba  near 
nodal  llnas .  This  paraonal  experience  was  aada  aora  objective 
In  a  study  by  Ihittll  (duttll  and  Oudaitla,  1966). 

NUttll  aaaaurad  tha  aaplltuda  of  tha  first  half -eye la 
of  tha  P  wave  at  numerous  stations  for  tha  Kodiak  Island 
aarthquaka  of  February  6,  196A.  These  aaplltude  data,  cor¬ 
rected  only  for  Instrumental  response,  ara  shown  In  Pleura  9.1. 
Tint  scatter  in  tha  data  U  considerable. 

Tha  PSult  plana  solution  of  this  particular  earthquake 
la  known  (Staudar  and  Bollinger,  1966).  The  theoretical  P 
wave  aaplltuda  radiated  by  a  double  couple  point  aourca  in 
a  homogeneous  medium  U  given  by  tha  relation 

•V  .0 

where  £  la  tha  density,  g  tha  P  wave  velocity,  R  distance 
along  tha  ray,  K(t)  tha  source  time  function,  and  a  and  y 
tha  space  coordinates  of  an  observing  point  with  respect  to 
the  mechanism  axes.  Ihittll  corrected  the  observed  amplitudes 
by  dividing  by  tha  factor  xy  appropriate  to  each  station. 


TM  rtaulu  art  ahoun  in  flfuf*  ^  .  Tha  «har»cur  of 
um  aynbola  in  tha  figure  la  rvlatod  to  tha  dUUnoi  of  • 
point  fron  •  nodal  11m,  or,  aqulvalantly,  to  too  angnltuda 
of  th«  factor  ay.  It  la  no  tod  that  tht  acattor  la  eonaldar- 
abl|  raducad,  with  um  axcaptlon  of  a  fra  polnta  for  Mich 
value a  of  ay  aro  loaa  than  0.1,  1  ,  polnta  nanr  nodal 

llnae.  By  perturbing  tho  nechanlaa  aolutlon  by  only  a  fan 
oagraaa  An  tha  dip  and  dip  dlreetlou  of  tha  nodal  planaa, 
a  aolutlon  wee  obtalnad  for  vhlch  a  van  thaoa  eaeoptlone 
dlaappaarad  and  for  which  tha  raalduala  in  obaarvad  vnraua 
oalculatad  anplltudaa  vara  a  nlnlnun. 

Tha  final  aolutlon  obtalnad  in  thla  way  la  ahown  in 
Figure  9  3-Tha  11m  dram  on  tha  figure  la  tha  aaplltude 
variation  expeoted  fron  gao— trie  apraadlng  only,  baaad  on 
tha  Joffray-Bullan  p  aava  trovol  tlna  eurvo  for  0.00B.  Tha 
aorraapondonaa  bat anon  tha  obaarvad  aaplltudea  aorroatod  for 
fooa*  Mohaalan  and  um  tnaoratlaal  anplltudaa  la  food.  Ona 
nay  aonoluda  that  tha  ahlaf  faetor  In  tha  aaattar  of  P  uava 


Thaoa  oboorvatlona  euggeat  tha  poaalblllty  of  aealrlng 
a  fooal  naohanlan  aolutlon  for  uhloh  aona  quantity,  any  tha 
vorlanao,  S,  of  tha  oboarvad  voroua  lanputad  anpUtudao  la 


Za  principle  this  1«,  of  course,  c  veil  MUbiUM  tech¬ 
nique,  thet  of  the  body-yve  squalls  tlon  procedure  which  use 
propoecd  bp  Pen  Wenshsm  et  cl.  (1966)  end  hoc  been  applied  to 
c  number  of  earthquakes  (Tong  end  Ben-dene  nee,  19661  Ben* 
R.nshoo  et  al.,  196&  Khsttrl,  1969)*  In  •  recent  paper 
Chandra  (1970)  bee  further  applied  the  method  of  equalise* 
tlon  of  P  wee  spectral  eapiltudee  to  determine  the  focal 
acehealao  of  flea  earthquakes,  end  has  compered  the  results 
with  solutions  from  the  first  notion  of  f  end  from  the 
polarisation  of  the  i  nave.  Be  found  close  agreement  in 
the  solutions  determined  Independently  by  the  three  pro* 

Za  the  application  here  proposed,  honeser,  there  la 
question  not  so  much  of  the  spectral  amplitudes  of  t  equal¬ 
ised  to  the  source  as  of  the  amplitude  of  the  first  half- 
eye  le  of  the  P  yea  ok e erred  in  the  time  domain.  That  is, 
the  aylitude  of  the  first  half-cycle  of  the  P  usee  is  con¬ 
sidered  to  represent  the  spectrum  of  the  yea  at  the 
predominant  period  of  the  P  yea.  This  is  a  stationary 
phase  apprealmetlon.  Tong  and  Ben  Bens  hem  first  used  the 
term  in  this  connection  and  applied  the  stationery  phase 
appronlamtlen  technique  to  the  SB  waves  of  the  Banda  lea 
earthquake  (tang  end  Ben  Bane  hem,  1966).  Chandra  (19T0b) 
has  applied  the  say  technique  systematically  to  the  P  yeas 
for  a  group  of  eight  earthquakes.  Bt  found  that  the  sta¬ 
tionary  phase  eolations  agreed  with  solutions  previously 
data— toad  by  equalisation  of  the  P  yve  spectra. 


A  stationary  phiit  •pprulMilon,  ■lnec  11  ultllMl 


only  the  first  nsif-cycl* .  or  cvtit  the  first  quarter  cycle 
(F  ism  onset  to  first  peek  or  trough)  sty  be  considered  to 
aeaple  the  direct  F  Arrival*  uncon teaina tod  Oy  converted 
us  vs  s  in  the  crust  benestA  im  ststions  or  In  Um  nolfAbor* 
hood  of  ths  source.  It  is  potentially  applicable,  there fore, 
to  she 1 low  focus  earthquakes  as  sell  as  to  the  isolated 
a  14ns la  of  F  waves  froai  deep  focus  shocks . 

Sven  after  equalisation  to  the  source,  F  wave  amplitudes 
will  very  In  absolute  value  depend ii*  upon  the  Magnitude  of 
the  earthquake.  In  order  to  cospare  equalised  aaplltudoa 
to  oslculstod  aaplltudea  a  noraalltlng  factor  suet  bo  lntro- 
duaod.  Par  this,  follouing  Jaroaoh  (1968),  «  lab 

A.,„  *  A  A.u* .  »  *•  (J) 


where  A^  ^  la  the  observed  asplltude  at  tbs  lth  station 

A1  calc  1*  M  oslculstod  asplltude  as  given  by  oq.(l) 
k  la  a  conatant  to  bo  do  to  mined 


Zn 

•i 


o4  la  the  error  in  the  1th  observation, 
order  to  nlnlnlaa  the  sun  of  the  aquaroo  of  the  a 
,  we  chooee  k  so  that 


k 


t 


(M 


A  value  of  k  nay  bo  so  obtained  for  any  orientation  of  the 
source  axes.  Put  a  search  pragma  la  latredeeed  In  order 
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to 


U»t  oriMUUw  of  U»  oouroo  ami  for  NtieU 


(5) 


is  o  minimum  A  boot  fitting  solution  is  thus  obtotnod  in 
o  in nor  quit*  onologouo  to  tho  procoduro  propoood  by  Bon* 
NmbIim  ot  •!.  (1964)  in  tho  •  quo  11  sot  ion  ond  oooporloon  of 
opootrol  u oplltudoo. 

An  llluotrotivo  oxonplo,  ond  o  prooonutlon  of  tho 
rooulto  in  0  plono  otnor  tnon  tno  googrophlo  plono,  iy  lip 
to  ounnorloo  botn  tno  itnod  ond  0  lloltotlon  ooonn  to  oil 
oquo  Union  toehnlquoo.  Plguro  4«  pro  00  n  to  tno  oonwntlonol 
P  woo  flrot  notion  ond  •  ww  polo  rl  tot  ion  ooononli  win* 
tloi  for  0  BnObo ldo  Mrtbquoto  of  Ootobor  n,  1969*  Xt  lo 
notod  tnot  tno  oslauthol  oowrogo  of  outlono  obout  tno 
oploontor  is  rotnor  good,  ond  tnot  ono  nodol  plono  lo  wll 
dotorn&nod  by  tno  P  woo  flrot  notion.  Tho  positions  of 
both  nodol  plows  aoy  bo  flood  by  0  boot  fit  to  tho  •  ww 
polorUtloo  dots. 

Xn  ordor  to  opply  tno  tqyoimtlon  of  tno  onplltudoo 
to  tho  foool  ophoro  by  tho  0  tot  ions  ry  phooo  oppro&lwtlon, 
tho  onplltodoo  ond  porlodo  of  tno  flrot  nolf-oyelo  of  tho 
P  ww  wro  do  ton  nod  ot  8)  outlono.  Thooo  oboorwtlow 
wro  tm  oquo Hood  to  tno  oouroo  by  using  Into  ooeount  tno 
offbot  of  tno  inotraoot  ond  of  U10  eniot.  ond  by  ooapon* 
Mtlng  for  goowtrlo  oproodlng  ond  oboorptlon  in  tno  nontlo. 
Tho  rooiiiu  of  tho  oooron  progn  oro  dloployod  In  Plguro  9.^. 


Tim  obMrvtd  •quallNd  MplUvtoi  an  proaontod  on  U»  loft, 
uw  cal  cult  tod  oopiitudoo  for  tut  boot  fitting  or  ion  tot  Ion 
of  tho  focol  oxoo  on  tno  rlgnt. 


Ao  hoo  boon  polntod  out  by  uooro  of  thla  technique, 
tin  grepho  of  Figure  *  >  oro  o  two  dlnonolonol  prooontotlon 
of  whet  in  roollty  lo  o  three  dlnonolonol  flguro i  tho 
irrogulorlt  loo  of  tbo  ooplltudo  dlogrono  oro  o  funotlon  not 
of  oolnuthol  dopondoneo  but  of  dlfforonoo  In  tno  tebe-off 
onglo  ot  tno  oouroo.  Per  exaaple,  Figure^  ohowo  tbo  vorlo- 
tlon  In  ooplltudo  oo  o  funotlon  of  oolnutn  oorroopondlng  to 
tno  grootoot  ond  to  tno  looot  toko -off  onglo  for  tno  par¬ 
ticular  oboonratlono  uood  In  tno  oxonl notion  of  tnio  earth- 
guano,  unile  oitnor  ounre  la  enooth,  it  la  roadlly  ooon 
that  if  royo  oorroopondlng  to  dlfforont  taho-off  ongloo  oro 
plot  tod  on  a  polar  graph  a  quite  lrrogular  flguro  oan  rooult. 
■moo  tho  otar-Uho  obopo  of  tho  oboarood  ant  oaloulatod 
ooplltudo  ploto. 

Mmo  for  ourfaoo  uavoo  It  la  noooooory  to  prooont  a 
radiation  pottom  la  a  googropnio  plana,  thla  dooa  not  ooon 
to  bo  noooooory,  or  ooon  of  particular  advantage  for  tho  re¬ 
duction  and  prooontotlon  of  body  wave  ooplltudo  data.  A 
praoontatlan  la  tamo  of  tno  fooal  a aaa  nay  bo  nco  Inattua 
tlwa. 

In  urao  of  opnarlool  ooordlnatoo  rofarrod  to  tbo 
oouroo  axoo,  tno  anpUtudo  of  tho  F  wave  on  tho  ourfaoo 
ophore  of  radluo  ■  (ooo  eq.  (l))  la  g&oon  by 


of  a 


»»«  1$  yr'/L  m  4 

4rr‘ * 


«6) 


vhtro  ^  It  neaaurad  in  lAt  i|  plant,  or  in  um  plant  of 
tha  force*,  ond  6  Is  araaured  from  in*  L  or  B  (null)  tala. 
TH»  variation  with  f  la  im  faalllar  four-lobad  F  wave 
radiation  pattam.  Both  f  and  0  an  readily  ra la  tad  to 
tna  sslauth  and  taxa-ofr  angles  of  a  given  rs;  through  tha 
direction  ooainaa  of  tna  avehanlaa  axaa  ralatlva  to  tha 
geographlo  axaa .  Tha  aaplltuda  alone  a  rap  asking  an  angle 
•  aoy  in  turn  ba  raduoad  to  tha  aaplltude  of  a  ivy  leaving 
a  aouroa  of  tha  aaaa  atrengl h  at  an  angle  t  •  7T/2  by 

divldint  tha  observed  aaplltude  by  air.2©,  in  thla  way  tha 
radiation  pattam  la  reduced  to  tha  xy  plana  or  plana  of 
tha  foraaa.  Or,  equl volant ly,  ona  views  tha  radiation 
pattam  froa  tha  B  axla,  looking  toward  tha  xy  plana. 

Figure  *7  presents  tha  anplltuda  data  of  tha  Hokkaido 
earthquake  in  thia  faahlon.  Ona  aaaa  iawed  lately  tha  a f foot 
of  tha  nodal  lina  and  tha  dagraa  of  closeness  of  fit  of  tha 
akaarvad  F  wave  anplltudaa  with  tha  calculated  aaplltudaa 
for  tha  Plantation  of  tha  focm  aa  lac  tad .  Ona  aaaa  aa  aa.l 
and  nara  graph! aally  than  in  othar  praaantationa,  tha  Halts 
tlon  in  tha  dlatrlbutlon  of  data  ovar  tha  focal  aphare. 

Ivon  with  good  geographic  dlatrlbutlon  or  obaarvatloao,  tha 
data  points  art  spread  over  barely  a  quarter  of  tha  xy  plana 


In  e«»«»  nowvir,  whoro  tin  distribution  in  tbo  x|  plono 


cover*  Hvtnl  qusdrsntft  (atrilo  slip  fruiting  in  •  vtrtlal 
plono)  •  prooontfttlon  ••  In  Plguro  M|  ftorvo  to  Idontlfp 
thft  plono  of  fruiting  and  tbo  diroot  ion  rupturo  propogotlom 
tbo  loboo  will  bo  ft*  yuoo  trie,  onlsrgod  in  tbo  diroot  ion  of 
tbo  rupturo. 

A  flnftl  rofrronoo  to  tbft  qyontlty  k  do  to  ml  nod  obooo 
(oq.  (4)).  if  uo  lot  X  bo  on  everogo  voiuo  of  tbo  oopiitudo 
on  tbo  foeol  opboro,  tbon 


b  nclog  1^ 


no  boot  N  M  log  (bXtftle) 
or  M  n  log  k  ♦  log  Xoa^c  . 

Tbo  ooiuo  of  I**l€  io  obtolaod  bp  tbo  roiotion 


Tbuo  tbo  uoluo  of  k  nop  bo  uood  to  dotoratao  o  nognitudo 
dorlvod  oinultonftouolp  fron  aonp  oboortotiono  of  tbo  P  nooo 
onplitudo  ond  oorrootod  for  tbo  offrot  of  tbo  itilitlun 
pottorn  of  tbo  oonroo. 

Tbo  ftbovo  technique*  «ro  bolng  oppllod  opt  tone  t  loo  lip  to 
tbo  fttudp  of  oortnquftKOft  In  tbo  Kurllo  Xftiondft  ond  In  Pom- 
leuftdor.  Tbo  rooultft  will  bo  bolpful  both  In  support  of  foeol 
ooebftnloo  studios  end  in  tbo  do  to  mine  t  Ion  of  doniasnt  noebonisno 
of  •  region.  nocoftsorp  for  efttobllobli*  roglonsl  noebsnlon  cor¬ 


rections. 
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Figure  2 


Figure  3 


Figure  4 


Figure  3 


Figure 


Figure  7 


Titles  to  Figures 

.  Amplitudes  of  first  half-cycle  of  the  P  wave, 
corrected  only  for  instrumental  response, 

Kodiak  Island  earthquake  of  February  6,  1964. 

.  Amplitudes  of  first  half-cycle  of  the  P  wa^e 
corrected  for  focal  mechanism  as  determined 
by  Stauder  (1956). 

.  Amplitude  of  first  half-cycle  of  the  P  wave 

corrected  for  adjusted  focal  mechanism.  Solid 
line  is  the  amplitude  expected  on  the  basis  of 
the  Jeffreys-Bullen  travel  time  tables. 

.  Conventional  first  motion  of  P  and  polarization 
of  S  mechanism  solution  for  the  Hokkaido  earth¬ 
quake  of  October  25,  1965. 

.  Observed  amplitudes  (left;  figure)  equalized  to 
the  source,  and  the  calculated  P  amplitudes 
(right  figure)  for  best  fitting  double  couple 
source,  earthquake  of  October  25#  1965. 

.  Variation  in  calculated  ?  wave  amplitude  for* 

take-off  ar.-'e-  _ ,.w  ...  .....  ^  _ 

and  least  of  ur  for  stations  represented 

in  Figure  5. 

.  Amplitude  data  for  the  Hokkaido  earthquake 
reduced  to  the  xy  plane . 
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10.  Source  ParawUn  of  larthguaites  from  Spectre 

of  R^rltlih  Mm 
A.  Udlaa 

A.  Ur^th  of  Mult  and  Velocity  of  Rupture. 

Ban-Menahe*  (1961)  haa  developed  the  expressions  for 
IM  far  flald  displacements  of  ftsylelgn  waves  radiating  fro* 
a  propagating  fault  of  flnita  dimensions.  The  aodal  of  the 
fault  la  a  aurface  of  length  b  and  width  d.  The  propagating 
fracture  la  repreaented  by  a  point  aource  propagating  along 
the  b  dimension  with  a  velocity  v  and  radiating  energy  aa 
It  propagatee.  Par  a  vertical  strike-slip  fault  and  a  prop¬ 
agating  double  couple  the  displacement  of  the  vertical  com¬ 
ponent  of  the  Rayleigh  wave  la  given  by 


The  effect  of  the  finite  dimensions  la  included  In  the 
factor  jCj|  (u»  »0).  To  laolate  this  factor  Ben-Menahe*  de¬ 
fined  the  directivity  function,  0,  aa  the  ratio  of  the  epec- 
trel  aeplltudea  of  the  aurface  wave*  correaponding  to  wave a 
leaving  the  aource  In  oppoalte  directions. 


ujttt  m  *«(*♦«). 


(3) 


Sine*  th*  only  factor*  which  dapand  on  6  In  oq.  (1)  in 
ain  XnAq  and  sin  0,  wa  can  fanarallaa  tha  axpraaalon  for 
tha  directivity  to  ba  tha  ratio  of  tha  apaotral  anplltudea 
at  two  a  tat  Iona  com  t  pond  lnq  to  ray  a  leaving  tha  aouroa 
with  an  arbitrary  ancle.  K  ,  batawn  thaw.  Thai' 

Minina  In  thla  function  occur  at  aaroa  of  tha  nunerator, 
that  la,  at  valuaa  of  tha  argunent  of  tha  alna  function 

■■■  (  •  (m%sj  a  (  #|a  2,  •••  • 

'  v  9 


Haxlna  ocou  at  saroa  of  tha  danonlnato r,  or  for 

*  wt  ,  ns  . . 


Tha  valuaa  of  tha  fraquaney  at  which  althar  aaxlna  or 
minim  occur  can  be  uaad  to  datarnlna  tha  fault  length  b. 
For  exaaple,  for  tha  kit  ortfar  axtraaa 


(5) 


The  dlmtlvlty  function,  D,  as  defined  originally 
(eq.  (l))  «i  uNd  vlili  R  or  0  mmi  of  eonMeuilvt  ordor 
roooriod  at  the  mm  station.  This  procedure  bIaIbImi 
tha  sources  of  error,  but  rodueoa  tha  uaa  of  the  aathod  to 
earthquates  sufficiently  large  to  generate  rat  urn  aurfaoa 
Mvat  or  Mvai  which  circle  tha  aar.n  more  than  orvea .  In 
theory  nothing  pravanta  using  dlraet  Rayleigh  or  Love  waves 
raoordad  at  different  stations  distributed  at  arbitrary 
but  differing  at  lout  ha  about  tha  ^Mcenter.  In  this  event 
tha  aora  general  directivity  function,  ,  (aq.  (a))  ap- 
pllaa.  Allowance  auat  ba  i  *da,  of  court a,  for  differences 
In  aploantral  dlatanea,  path  of  propagation,  ate  This 
prooadura  nates  It  post  lb  la  to  apply  tha  aathod  to  earth¬ 
quakes  ana  liar  than  M  •  7. 

Tha  advantage  of  applying  tha  aora  generalised  expres¬ 
sion  for  tha  dlreetlvlty  la  that  stations  can  ba  aeleeted 
at  which  tha  Rayleigh  waves  are  wall  raoordad  without  tha 
constraint  of  being  100°  apart  in  aslauth.  Since  tha 
directivity  function  it  noraa Used,  tha  amplitudes  in  tha 
apaotral  ratios  auat  ba  reduced  to  a  coaaon  dlatanea.  This 
la  dona  through  tha  aquation 

A,  •  At(eKf 

Ae  an  axaapla,  tha  aathod  hare  proposed  la  applied  to 
tha  Aleutian  Island  earthquake  of  July  A,  1966.  Tha  two 
nodal  planes  of  f  are  assumed  to  ba  known,  and  are  taken 


froa  dtau tier  (**'*rj.  «c .  .  1  anows  Ww  orlantetlona, 

aa*  and  hb' .  of  w.«*c  p.*:««*  and  also  indicate*  the  a  tat  Iona 
and  UMlr  dlatritution  about  tn*  aourca  at  which  tha  direct 
Rayleigh  w«v«a  ware  uaed  in  tha  anal/ala .  The  Rayleigh 
wave  a  at  thaaa  a  tat  Iona  ware  digitised  and  a  Courier  ana  lye  la 
■ada  of  tha  digitised  data.  In  thia  procaaa  care  auat  ba 
taken  to  aalact  a  suitable  tine  window,  for  it  la  Important 
to  exclude  froe  tha  analysis  an/  interference  froa  later 
arrive  la.  lateral  refractions,  ate. 

Oiven  tha  afarctra,  aa  aeny  combinations  aa  possible  of 
apeetral  ratio*  are  foiaed.  By  application  of  equation  ()), 
repeated  froe  above 

u _ 1*  »  Sa -  (5) 

JJT  ^  -  «•*(•♦«}]  A_(  ^  *<«*♦) 

a  value  of  b  can  be  obtained  fro*  the  first  nasleun  or  rirst 
alnleue  of  any  one  of  thaaa  apeetral  ratios.  Zn  so  doing 
It  la  necessary  to  aaauae  a  value  of  v,  the  velocity  of 
rupture.  In  fact,  when  ualng  only  the  amplitude  portion  of 
the  apactra  there  la  no  Independent  way  of  date  raining  the 
fracture  velocity,  v.  However,  if  narvy  ratioe  are  uaed  for 
the  calculation  of  b,  we  nay  aaauae  different  valuea  of  v 
and  froo  thaaa  aalact  that  v  which  glvea  the  leaat  atandard 
error  for  tha  naan  value  of  b.  Zn  thia  wanner  we  obtain  a 
alaultaneoua  aatlaate  of  both  b  and  v. 


177 


A  by-product  of  the  set  hod  when  %  sufficient  number  of 
stations  s round  the  epicenter  are  used  is  the  determination, 
through  the  proper  choice  of  0  of  the  plane  of  faulting  from 
the  two  nodal  planes  and  the  direction  of  fracture. 

Pigs. 10. d  k  It.;  show  two  examples  of  the  spectral 
ratios  for  the  earthquake  of  July  *♦,  1^66.  The  points  in 
the  figures  are  observational  values  of  the  spectral  ratio 
for  the  pairs  of  stations  indicated.  The  solid  lines  are 
the  computed  directivity  function  for  the  fault  length, 
fracture  velocity,  and  direction  of  fracture  assumed  for 
the  beat  fit  to  the  data.  While  mere  Is  some  variation  in 
the  observed  from  the  computed  spectral  ratio,  the  minima 
and  max lam  are  well  defined.  For  this  Aleutian  Island  earth¬ 
quake  the  NS  striking  plane  is  selected  as  the  plane  of  fault¬ 
ing,  with  the  fracture  propagating  southward.  This  is  in 
keeping  with  the  Interpretation  of  mis  ear cnquake  as  an 
arc-arc  transform  fault. 

Pig. 10. A  shows  a  similar  application  to  an  earthquake 
in  tha  Azores.  The  focal  mechanism  was  determined  by  Sykes 
(1967).  The  earthquake  la  atrlae-sllp,  with  epicenter  on 
a  fracture  zone. 

b.  Seismic  Moment. 

The  eelsmlc  moment  la  calculated  following  the  work  of 
Akl  (1966)  end  the  theoretical  results  of  Ben-Kenahem  and 
Ilarkrlder  (1964)  from  the  spectral  amplitudes  of  the  vertical 
component  of  the  Rayleigh  waves  at  periods  sufficiently  large 
that  tha  point  source  approximation  of  the  source  is  valid. 


17” 

Assuming  a  step  function  fcr  the  source  time  function  the 
value  of  the  moment  la  given  by 

*"  NrJXW| 

U(vm)  la  the  apectral  amplitude,  r  the  distance,  ^  the  atten¬ 
uation  coefficient,  CR  the  phaae  velocity,  Nrz  the  Rayleigh 
wave  singlet  transfer  function  and  %(Q)  the  radiation 
pattern  function  as  defined  by  Ben-Menahem  and  Harkrlder. 

For  the  earthquakes  studied  here  we  calculated  the  moment 
at  two  values  of  the  period,  100  and  50  seconds.  For  these 
periods  the  uncertainties  in  the  assumed  phase  velocity  and 
attenuation  coefficient  are  small.  The  results  are  given 
in  Table  1. 


C.  Stress  Drop,  Seismic  Energy,  and  Average  Dislocation. 

Values  for  the  stress  drop,  ,  and  average  dislocation, 
w  ,  have  been  found  from  the  above  determined  values  of  the 
moment  and  fault  length  using  the  relations 


—  -  M» 

A  "  TTbT 


1 


where 


(8) 

(9) 

(10) 


The  value  of  was  taken  as  3.3  x  10"  and  the  fault  width 
equal  to  10  km  for  the  Aleutian  earthquake  and  5  km  for  that 


in  the  Azores.  The  results,  again,  are  given  in  Table  1. 
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Figure  1 


Figure  2 


Figure  3 

Figure  4 


Titles  to  Figures 

Orientation  of  nodal  planes,  AA 1 ,  33',  and 
azimuthal  relation  of  stations  selected  for 
surface  wave  analysis,  Aleutian  Island  earth¬ 
quake  of  July  4,  lSoo. 

Spectral  ratio  for  station  pair  ATU/AFI. 

Plotted  points  correspond  to  observational 
value:  solid  line  to  the  theoretical  spectral 

ratio  for  the  values  of  b  and  v  indicated, 
Aleutian  Island  earthquake  of  July  4,  1966. 

Spectral  ration  for  station  pair  ANP/BHP, 
Aleutian  Island  earthquake  of  July  4,  1966. 

Spectral  ratio  for  station  pair  TRI/IST, 

Azores  earthquake  of  July  4,  1966. 
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11.  Southeast  Ml as our 1  Seismic  Network 
Ram  Agrawal 

The  Immediate  objective  of  this  research  has  been  to 
study  the  selamlclty  of  the  New  Madrid  zone  on  the  basis  of 
small  earthquakes  (2.4  <  <3-5)  recorded  by  portable  short- 

period  sensitive  seismographs.  The  long-range  goal  is  to 
determine  the  relationship  of  the  spectral  properties  of  the 
seismic  waves  to  the  earthquake  parameters,  especially  fo?al 
depth  and  magnitude. 

Oeology  of  the  New  Madrid  Area. 

The  area,  geologically, Is  essentially  confined  within 
the  borders  of  the  northern  part  of  the  Mississippi  embay- 
ment.  This  embayment  Is  a  broad  northward  extension  of  the 
Oulf  Coastal  plain.  This  area  Is  bounded  on  the  west  and 
north  by  the  Ozark  Uplift,  and  on  the  east  by  the  Illinois 
Basin.  Both  of  these  major  features  are  made  up  of  sedi¬ 
mentary  rocks  of  Paleozoic  age.  Structurally,  the  embayment 
Is  a  large  downwarped  region  of  Paleozoic  rocks  opening  to 
the  Oulf  of  Mexico,  filled  to  the  level  of  present  surface 
with  deposits,  mostly  unconsolidated,  ranging  In  age  from 
Cretaceous  to  recent.  A  schematic  geological  cross-section 
W-E  at  latitude  3o  N  Is  shown  here. 
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Location  of  Stations 

A  three-station  seismic  network  has  been  operating  In 
western  Tennessee  since  October  1969,  with  the  following 
specifications: 

Samburg(S) 

Long.:  89.303°W 
Lat.:  36.4l6°N 
Elevation:  138  meters. 

The  station  Is  on  Tertiary  rock.  Recording  Is  on  mag¬ 
netic  tape,  played  back  In  the  laboratory  on  a  Brush  oscillo¬ 
graph.  Three  Instruments;  one  vertical  component,  one  hori¬ 
zontal  N-S  and  another  horizontal  E-W,  are  operating.  The 
natural  period,  and  peak  magnification  of  the  instruments 


are  as  follows: 
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Instrument 

Natural 

Period 

Amplifier- 

Gain 

Peak 

Magnification  at 

Frequency 
in  cps 

Vertical 

1  sec 

150 

33K 

13.5 

500 

106K 

13.5 

Horizontal 

(N-S) 

1  sec . 

150 

63K 

11.0 

500 

208K 

11.0 

Horizontal 

(E-W) 

1  sec 

150 

47K 

12.0 

500 

146K 

12.0 

The  noise  level  at  Samburg  has  been  found  to  be  5  millimicron 
peak  to  peak  at  11.0  cps. 

Lassiter  (L) 

Long.:  89.352°W 

Lat.:  36.330°N 

Elevation:  108  meters. 

This  station,  about  11  km  distant  from  Samburg  on  Ter¬ 
tiary  rock,  consists  of  one  vertical  component  seismograph 
with  natural  period  0.4  sec  and  peak  magnification  is  110K 
at  2.5  cps.  Recording  is  on  an  Autocorder. 

Tlpton^llle  (T) 

Long.:  89 . 507°W 
Lat.:  36.443°N 

Elevation:  85  meters . 

A  vertical  component  seismometer  with  natural  period 
1  sec  is  located  on  alluvium  at  a  distance  of  18  km  from 
Samburg.  Recording  is  on  magnetic  tape. 

In  addition,  another  station  is  operated  in  Missouri 
about  80  km  from  the  three  in  Tennessee: 


Greenville  (G) 

Long.:  90.  }9*>°W 

Lat.:  J7.0‘);J°N 

Elevation:  loo  meter 

A  vertical  component  seismometer  with  period  1  sec  and 
peak  magnification  10^  at  2.o  cps  recorded  on  an  Autocorder 
is  operated  there. 

Epicenter  Location  Technique 

Prom  a  model  in  which  all  stations  are  assumed  to  lie  on 
a  horizontal  plane  and  all  crustal  layers,  of  constant  veloc¬ 
ity.  are  assumed  to  be  bounded  by  horizontal  planes,  two 
computer  programs  have  been  developed  for  determining  hypo- 
centers: 

(1)  Using  Pg-arrlval  times  -  t^,  t^,  t^,  t^  at  the  four 
stations  (Westphal  &  Lange  1967).  Where  tj  -  time 
of  arrival  of  Pg  at  the  nearest  station. 

(2)  Using  Pg-arrlval  times  at  three  stations  and 
(Sg-Pg)  Interval  at  the  nearest  station. 

Both  programs  compute  focal  depth,  but  only  in  a  few 
cases  were  meaningful  depths  obtained.  Seismic  events  re¬ 
corded  from  March  to  May  1970  have  been  studied. 

In  Table  1-A  and  B  epicentral  coordinates,  origin  time 
and  focal  depths  are  given. 

The  following  velocity  structure  (McEviily,  1964)  was 


used  to  locate  the  epicenters: 
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VP 

vs 

Density 

Thickness 

6.1  km/sec 

3.5  km/sec 

2.7  gm/cc 

11  km 

6.4  km/sec 

3.68  km/sec 

2.9  gm/cc 

9  km 

6.7  km/sec 

3.94  km/sec 

2.9  gm/cc 

18  km 

8.15  km/sec 

4.75  km/sec 

3.3  gm/cc 

epicenters  were 

located  on  the 

basis  of  one 

of  the  follow 

ing  alternate  assumptions: 

1.  depth  less  than  11  km 

2.  depth  between  11  km  and  20  km. 

The  results  from  the  first  assumption  are  given  in  Table  1A 
and  those  from  the  second  in  Table  TB. 

Interpretation  of  the  Epicenter  Plot. 

Epicenters  thus  obtained  were  plotted.  Most  of  the  epi¬ 
centers  lie  from  35°  to  37°N  and  from  89°  to  90°W.  There  are 
two  clusters  of  epicenters;  one  near  the  New  Madrid  bend  of 
the  Mississippi  River  and  another  at  the  border  of  Missouri 

and  Arkansas.  However,  both  clusters  of  epicenters  lie  along 

o  o 

a  line  running  from  S  35  W  to  N35  E,  which  is  more  or  less 
parallel  to  the  course  of  the  Mississippi  River  as  shown  in 
Pig.  11.1.  The  epicenters  in  Pig.  11.1  are  under  assumption 
1.  above.  This  line  of  epicenters  is  in  good  agreement  with 
the  line  of  epicenters  based  on  historical  data  for  1811-1964 
(Heinrich,  1941,  Stauder,  1965*  McGinnis,  1963). 

The  above  line  of  epicenters  is  also  in  good  agreement 
with  the  direction  of  faults  located  by  Mateker  (1968)  using 
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gravity  and  well-log  data.  The  events  do  not  fall  on  any  of 
Mateker's  faults,  but  in  between  the  two  that  bound  the 
upper  end  of  the  embayment. 

Depth  of  Focus 

For  two  events,  depth  of  the  source  has  been  found  to 
be  around  12  km  and  that  for  one  event  17  km.  These  are 
expected  values  in  this  region.  In  three  cases  tentative 
values  of  focal  depth  of  100-200  km  were  found.  This  is  an 
unexpected  result  and  must  be  verified  by  further  analysis . 
The  tentative  spatial  distribution  of  these  six  events  along 
the  line  of  epicenters  (S35°W  to  N35°E)  is  shown  in  Fig. 
11.2. 

Any  definite  relationship  between  the  shallow  and  deep 
focus  events  cannot  be  given  because  of  the  limited  data. 

For  most  of  the  events,  computed  depth  of  the  focus  has  been 
found  to  be  Imaginary.  Thus,  effort  is  still  going  on  to 
develop  a  method  to  compute  the  focal  depth  of  such  events . 
It  has  been  Judged  that  only  if  the  epicenter  is  within  the 
network  and  also  very  close  to  one  of  the  stations,  the 
focal  depth  computed  is  reliable. 

Time  Residual  in  Pg-Arrival. 

Residual  in  Pg-arrival  times  using  three-station  tech¬ 
nique  of  epicenter  location  has  been  found  to  be  about  + 

0.3  sec.  But  in  four-station  technique  it  is  about  +2.0 
secs,  because  of  the  limited  control  at  Greenville  station, 
which  is  about  80  km  from  the  tripartite  network. 
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Magnitude  Determination. 

There  was  no  calibration  function  available  to  calculate 
the  magnitude  of  local  earthquakes  in  the  New  Madrid  area. 

A  relation  for  the  calculation  of  magnitude  has  been  developed 
using  Pg. 

Development  of  the  Magnitude  Calibration  Function. 

Let  us  assume  the  magnitude  relationship  as  follows: 
mb  *  a  +  log  £  (-g-fg— )  +  b  log  D  (km)  (1) 

where 

A  -  ground  motion  in  microns 
T  -  period  of  Pg 

D  -  hypocentral  distance  in  km 

Ip  o 

D  -  hd  +  £>  d 

h  -  focal  depth 

epicentral  distance  in  km. 
a  and  b  are  constants 

a  -  varies  from  event  to  event 

b  -  depends  upon  the  local  geological  conditions. 

A 

Thus,  for  about  sixteen  earthquakes  log  ip  and  log  D  are 
calculated  at  Samburg  and  Greenville,  assuming  focal  depth  to 
be  12  km,  and  are  shown  in  Table  II. 

A 

Then  log  p  vs  log  D  curve  has  been  plotted.  The  shapes 
of  the  straight  lines  thus  drawn  for  each  event  are  measured 
as  shown  in  Table  II.  The  average  slope  turned  out  to  be  tan 
65°;  i.e.  b  =  tan  63°  *  2.15. 

h  m 

The  magnitude  of  the  event  of  March  27,  1970,  03  44  , 

was  calculated  to  be  3*3  by  Dr.  0.  Nuttli  using  Pn.  Taking 
this  earthquake  as  a  standard  one,  a  line  of  slope  of  2.15 


< 
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was  drawn  passing  through  the  point  (log  ^  ,  log  D)  for 

A 

Samburg  as  shown  in  Pig.  11.3.  This  line  cuts  log  axis 
at  I.63  where  log  D  =  1.0.  Thus  on  substitution  the  values 
of  mb,  log  fjr  and  log  D  in  equation  (1), 

3.5  =  a  +  1.63  +  2.15  x  1.0 
a  =  -3-78  +  3.50  =  -0.28. 

Therefore,  magnitude  relationship  on  the  basis  of  16  earth¬ 
quakes  using  Pg  becomes 

mb  =  -0.26  +  log  £  ("sec0-')  +  2.15  log  D  . 

Magnitude  of  all  the  events  were  calculated  using  this  rela¬ 
tionship  and  are  shown  in  Table  II.  The  magnitude  varies  from 
2.4  to  3.5. 

Justification  of  the  Magnitude  Formula. 

MM 

The  magnitude  of  the  event  March  27,  1970,  02n  21  ,  was 
calculated  to  be  2.5  by  Dr.  0.  Nuttli,  and  by  this  relation¬ 
ship  also  it  has  been  found  to  be  2.5.  Hence  the  above  mag¬ 
nitude  relationship  seems  to  be  reasonable;  however,  for 
better  magnitude  relationship  more  data  are  needed. 

Discrimination  of  Earthquake  and  Blast. 

At  about  100  km  from  the  tripartite  network  there  are 
mining  districts  where  frequent  blasting  is  done.  Thus  the 
discrimination  of  local  earthquake  and  blast  in  the  New  Madrid 
area  is  very  difficult.  However,  the  following  criteria  have 
been  used  for  discrimination: 

1.  On  the  basis  of  epicenter  location.  Epicenters  of 
some  of  the  events  have  been  found  to  be  in  the  known 
mining  districts;  for  example,  the  events  of  Apr.  19, 
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1970,  10h  ll”*  and  Mar.  28,  1970,  07h  38m  . 

2.  On  the  basis  of  appearance  of  seismograms .  The  seismo¬ 
grams  of  strip  mine  explosions  and  local  earthquakes 
have  been  shown  in  Pig.  11.4.  As  a  general  rule,  rec¬ 
ords  from  industrial  explosions  in  the  Illinois -Missouri 
area  feature  a  well  developed,  short-period  surface  wave- 
train,  such  as  that  seen  in  Pig.  11,4  b,  which  may 
easily  be  recorded  to  distances  of  200  km  and  greater. 

The  records  from  small  earthquakes,  on  the  other  hand, 
are  notably  lacking  in  the  surface  wave  trains. 

3.  On  the  basis  of  ratio  of  Pg  to  Sg.  Ratios  of  Pg  to  Sg 
ha\e  been  calculated  for  all  the  events  and  are  shown 
in  Table  II.  It  has  been  found  that  a)  Pg/Sg  ^  1  for 

earthquake,  b)  Pg/Sg  >  1  for  blast.  For  example, 

,  _  h  m 

Pg/Sg  -  1.71  for  the  event  of  Mar.  28,  1970,  07  38  . 

Seismicity  of  New  Madrid  Zone. 

Prom  the  recording  of  Mar. -May,  1970,  it  has  been  found 
that  in  a  month  six  events  occurred.  In  the  following 
table  the  number  of  events  of  a  particular  magnitude 


are  shown: 

No.  of  events  in  three  mo.  mb 

2  2.4 

3  2.5 

3  2.8 

1  3.0 

1  3.2 

2  3.5 
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The  range  of  magnitudes  in  the  present  data  set  is 
too  small  to  permit  the  determination  of  a  meaningful  slope 
for  the  log  N  vs  mb  curve. 

Conclusion. 

On  the  basis  of  the  earthquakes  occurring  from  1962-69 
as  reported  by  USCGS,  it  has  been  concluded  that  there  are 
only  five  events  per  year  in  the  Southeast  Missouri  area. 
However,  this  network  recorded  about  six  events  per  month 
during  March-May,  1970. 

In  view  of  this  fact,  though  at  present  the  results  are 
not  complete,  it  is  anticipated  that  continued  recording  will 
result  in  a  more  accurate  picture  of  seismicity. 
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Table  1-A 


April  08 
May  08 


22 

23 

28 

30 


Origin  Time  Eplcentral  Coordinates  Focal  Depth 

H  M  S  Long .  in  Deg .  W  lit,  in  Deg.  N  in  Km  Remark 


22  42  31-39 

90.13 

35.93 

X 

* 

03  44  29.08 

89.48 

36.55 

11.88 

* 

07  38  30.31 

98.32 

28.46 

X 

* 

06  03  50.26 

89*53 

36.53 

X 

* 

03  41  44.54 

90.01 

35.67 

81.47 

* 

08  59  46.94 

89.89 

35.98 

X 

* 

10  29  48.63 

88.43 

35.92 

X 

* 

03  52  33.79 

90.02 

36.06 

X 

* 

23  50  37.69 

89.63 

36.16 

11.22 

* 

08  04  16.29 

90.06 

35.48 

X 

* 

11  24  03.97 

89.82 

36.63 

X 

* 

23  39  10.83 

89.55 

36.43 

5.22 

* 

05  02  58.75 

89.05 

37.04 

150.86 

* 

March  27 

02  21  34.51 

89.59 

36.58 

X 

# 

03  44  27.79 

89.51 

36.60 

18.10 

# 

28 

07  38  29.47 

95.82 

30.77 

X 

# 

April  11 

00  59  55.62 

89.61 

36.26 

X 

# 

14 

04  39  10.76 

89.44 

36.33 

11.83 

# 

19 

10  11  23.37 

86.77 

40.11 

X 

# 

Vp  -  6.1  tan/sec  X  -  Depth  Imaginary 

Vs  »  3.5  km/sec  *  -  Three  Station  Data  Technique 

#  -  Four  Station  Data  Technique 
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Table  IB 


Date  1970 

Origin  Time 

H  M  S 

Epicentral 
Long,  in  Deg. 

Coordinates 
y  EaF." "fri  'Peg' ."  "N 

Focal  Depth 
in  km 

Remark 

March  18 

22  42  31.38 

90.21 

35.93 

X 

* 

28 

02  21  34.91 

89.58 

36.58 

X 

* 

03  44  29.07 

89.49 

36.57 

11.44 

* 

28 

07  38  30.29 

96.67 

29.98 

X 

« 

April  08 

06  03  58.99 

89.64 

36.35 

184.35 

* 

11 

00  59  55.54 

89.40 

36.65 

X 

* 

14 

04  39  13.44 

89.41 

36.42 

X 

* 

May  08 

03  41  58.44 

90.02 

35.57 

81.08 

* 

08  58  46.89 

89.89 

35.92 

X 

* 

10  29  48.62 

88.28 

35.85 

X 

* 

22 

03  52  33.76 

89.16 

35.67 

X 

« 

23  50  37.65 

89.6O 

36.11 

13.03 

* 

23 

08  04  16.28 

90.08 

35.35 

X 

* 

11  23  58.84 

90.39 

35.84 

X 

* 

28 

23  39  10.71 

89.51 

36.41 

17.12 

* 

30 

05  02  58.61 

89-95 

37.09 

157.18 

* 

March  27 

02  21  35.05 

89.58 

36.57 

X 

# 

03  44  28.48 

89.50 

36.59 

14.34 

# 

28 

07  38  30.13 

96.10 

30.49 

X 

# 

April  11 

00  59  56.47 

89.59 

36.28 

X 

# 

14 

04  39  12.01 

89.43 

36.36 

X 

# 

19 

10  11  23.29 

86.50 

40.40 

X 

Origin  Time  Frequency  Samburg  Greenville  b  =  tanG 

Date  1970  H  M  S  S-P  Pg/Sg  in  cps  log  A/T  log  D  Log  A/T  log  D  6  in  deg. 
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Figure  11.4.  Comparison  of  records  from  BILLIKEN  and 
other  small  seismic  sources  of  the  Missouri  area, 
a}  (top)  BILLIKEN  record  at  DeSoto,  Missouri, 

b)  (middle)  strip  mine  explosion,  also  recorded 
at  DeSoto. 

c)  (bottom)  small  local  earthquake  recorded  near 
Poplar  Bluff. 
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